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ABSTRACT 


Standard operating procedure at the Kennedy Space Center (KSC) for the 
Space Shuttle Program requires the storage and transfer of substantial quantities of 
liquid hydrogen (LH 2 .). Vaporized liquid, routinely lost during these transfer 
operations, is vented to the atmosphere or burned in the bum pond, and represents a 
significant fraction of the total hydrogen-fuel used for each launch. This report 
describes a procedure which uses metal hydrides to capture some of this low 
pressure (<1 psig) hydrogen for subsequent reliquefaction. Of the five normally 
occurring sources of boil-off vapor the stream associated with the off-loading of 
liquid tankers during dewar refill was identified as the most cost effective and 
readily recoverable. The design, fabrication and testing of a proof-of-concept 
capture device, operating at a rate that is commensurate with the evolution of 
vapor by the target stream, is described. Liberation of the captured hydrogen gas at 
pressures >15 psig at normal temperatures (typical liquefier compressor suction 
pressure) are also demonstrated. A payback time of less than three years is 
projected. 


I 

I 


I. EXECUTIVE SUMMARY 


A. Purpose 


The volume of liquid hydrogen which will be used during the late 1980’s 
by the Space Shuttle program at the Kennedy Space Center is expected to exceed 10 
million gallons per year. Of this 10 million gallons only about 7 million will actually 
be used as propellant. The balance, over 3 million gallons per year, will be lost as 
boiloff in various transfer and storage operations at launch pads 39A and B. 

The highest valued use for hydrogen boiloff identified during the 
preparation of the unsolicited proposal which led to this work was as feed stock for 
reliquefiers. Unfortunately, the boiloff rates during the intermittent loss periods 
are too high for economically sized equipment to reliquefy in "real time". The loss 
rates are also too high for conventional gas compressors and storage equipment to 
capture the boiloff vapor. 

Metal hydride hydrogen absorbers are fast enough to capture the boiloff 
hydrogen during many of the loss episodes at LC-39. Hydrides are fairly expensive 
however, so the only losses which may be economically captured are those which are 
large and occur frequently. This leads to many cycles of use for the hydride and a 
large credit for hydrogen capture. 

The most attractive use for hydrides identified in this program is for 
capturing the losses which occur during the transfer of liquid hydrogen from truck 
and rail cars. These "off-loading" operations will occur about 170 times per yeard). 
Other opportunities to capture boiloff are apparent which bring the total annual 
cycle number to about 200. Approximately 1.5 million gallons can be recovered 
from this source which will result in an annual savings >$2 million (assuming 
$ 1.50/gallon). At this rate of use, a metal hydride absorber will pay for itself in a 
few years of use (2-3 years depending on reliquefaction costs). 

B. Approach 

In the course of this program a "proof-of-concept" unit utilizing the 
unique characteristics of metal hydrides was designed, fabricated and tested. The 
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device was used with low pressure (atmospheric) hydrogen gas to demonstrate the 
technical feasibility of recovering hydrogen normally vented at the KSC. 

C. Results 


The test results indicate that a metal hydride capture/storage system: 

. Can absorb low pressure (<2.5 psig) hydrogen at a rate that is 
commensurate with the rapid boiloff of hydrogen vapor 
experienced during the LH 2 tanker "off-loading” operation. 

. Will interface with existing equipment at LC-39 with minimum 
impact on present operating procedures and be compatible with 
KSC safety practices. 

. Can release the stored hydrogen at normal temperatures and at 
pressured >2 atma for reliquefaction at a rate that is synchronized 
with the LH 2 delivery schedule. 

. Provides sufficient economic incentive to warrant consideration 
for further development of the concept at the KSC. 

The tests confirmed the validity of the initial design configuration. A 
full-scale storage system (to absorb one-hour boiloff) will require 290,000 lbs of 
metal hydride. The hydride will be contained in 1-1/8 IN O.D. copper tubes with a 
flexible filter (gas distributor) on the axial center line for the full length of the coil. 
The tubes are spirally coiled with a minimum diameter of 2 feet and a maximum 
diameter of 10 feet. 73 coil layers are stacked in each of the five 10 ft. high 
vessels. See Tables 1-3 for three of the size options considered. Table 4 is a listing 
of the program that was used to generate the weight and size requirements of the 
full-scale system, and Table 5 is a compilation of thirty-six computer runs that 
represent the range of reasonable sizes. 
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INPUT DATA 


TABLE Is 


TABLE 2: 


TUBE DIAMETER • 1.123 IN 
COIL DIAMETER — HA* I«UH - 10 FT 
SPACE BETWEEN COILS - .3 IN 
HYDRIDE CAPACITY - 1 X 
HA* I hum STACK HEIGHT * 15 FT 


NUMBER OF COILS - 29.34 
HYDRIDE WEIGHT - 764.16 LBS 


WALL THICKNESS - . 03 IN 

HINIHUM - 2 FT 

BETWEEN LAVERS • .3 IN 
PACKING DENSITY - 60 % 

AXIAL FILTER OD .223 IN 

layer totals •••••• 

COIL LENGTH - 336.3 FT 
HYDRIDE VOLUME - 1.824 FT3 


SYSTEM TOTALS 


COIL LAYERS PER STACK - 110 NUMBER OF STACKS REQUIRED -4.04 

TUBING LENGTH - 247603 FT INTERNAL VOLUME - 1332.6 FT3 


HYDRIDE VOLUME - ail. 36 FT3 22.99 M3 

HYDRIDE WEIGHT • 340000 LBS 134343.43 KGS 

HYDROGEN STORED - 3400 LB9 1343.43 kQB 

HYDROGEN STORED * 618181.82 SCF l?3i2.23 SCH 

HYDRIDE COST • 63090909 SYSTEM COST • 64327272 

PAYBACK • 61.8MM/YR - 2.4 YEARS 


Computer Run for 10 ft. diameter coil, 15 ft. stack height and 1% 
hydride storage capacity. 


INPUT DATA 


TUBE DIAMETER - 1.123 IN 
COIL DIAMETER — HA* I MUM - 
SPACE BETWEEN COILS - .3 IN 
HYDRIDE CAPACITY - 1.2 1 
MAXIMUM STACK HEIGHT - 10 FT 


WALL THICKNESS - .03 IN 

0 FT MINIMUM * 2 FT 

BETWEEN LAYERS - .3 IN 
PACKING DENSITY - 60 X 
AXIAL FILTER 00 .223 IN 

LAYER TOTALS ••*••• 


NUMBER OF COILS - 29.34 COIL LENGTH - 336.3 FT 

HYDRIDE WEIGHT - 764.16 LBS HYDRIDE VOLUME - l . S24 FT3 


*•*••• SYSTEM TOTALS 


COIL LAVERS PER STACK - 73 
TUBING LENGTH - 206337 FT 


NUMBER OF STACKS REQUIRED -3.08 
INTERNAL VOLUME - 1127.16 FT3 


HYDRIDE VOLUME - 676.3 FT3 

HYDRIDE WEIGHT - 293333.333 LBS 

HYDROGEN STORED - 3400 LBS 

HYOROGCN STORED - 610181.02 SCF 
HYDRIDE COST - *2373737 
PAYBACK • 61.BHM/YR • 2 YEARS 


19. 14 «3 

12B7B7.es KGS 

1343.43 KGS 

17312.23 SCH 

SYSTEM COST • *3*06O«0 


Computer Run for 10 ft. diameter coil, 10 ft. stack height and 1.2% 
hydride storage capacity. 


- 3 - 


INPUT DATA 


TU®€ DIAMETER • 1.123 IN WALL THICKNESS ■ .03 IN 

COIL DIAMETER — MAXIMUM - 12 FT MINIMUM - 2 FT 

SPACE BETWEEN COILS - .3 IN BETWEEN LAYERS - .3 IN 

HYDRIDE CAPACITY - 1.4 X PACKING DENSITY - 60 *4 

MAXIMUM STACK HE I GHT - 13 FT AXIAL FILTER OD .223 IN 

****•♦ LAYER TOTALS •••••• 


NUMBER OF COILS - 36. 92 COIL LENGTH • 811.37 FT 

HfDRIOE WEIGHT • U14.A uBS HYDRIDE VOLUME - 2. 636 FT3 

****** SYSTEM TOTALS •**••• 

COIL LAYERS PER STACK * 110 NUMBER OF STACKS REQUIRED *1.98 

TUBING LENGTH - 176862 FT INTERNAL VOLUME - 963.41 FT3 

HYDRIDE VOLUME - 379.23 FT3 16.41 m3 

HYDRIDE WEIGHT - 24283?.143 LBS 110389.61 KGS 

HYDROGEN STORED - 3400 LBS 1543.43 KGS 

HYDROGEN STORED - 618181.82 SCF 17312.23 SCM 

HYDRIDE COST - 92207792 SYSTEM COST • 93090909 

PAYBACK 8 91, 8MM/ YR • 1.72 YEARS 


TABLE 3: Computer Run for 12 ft. diameter coil, 15 ft. stack height and 
1.4% hydride storage capacity. 

D. Recommendations 


The next phase should include the following tasks: 

1) The device fabricated in the performance of this contract should be 
tested using hydrogen from LH 2 boiloff. 

2) An engineering analysis should be made of the various options available 
for hydrogen capture and storage for the purpose of identifying the most 
cost effective approach. 

3) A single, full-scale component should be fabricated to operate in-situ at 
LC-39 to field test LH 2 boiloff capture capability. 

4) Based on the results of all previous work, a full-scale system should be 
designed. 

5) An RFQ should be issued to solicit bids for construction of the full-scale 
capture system. 
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1 HOME 

10 D« * CHR* (4): I* - CHR* <9) 

15 REM 

20 REM **♦* THIS IS COIL STACK. **** 

25 REM 

30 REM THIS PROGRAM CALCULATES THE COIL LENGTH, VOLUME AND HYDROGEN CAPACITY F 
OR A SPIRALLY WRAPPED LAYERED STACK 

40 REM INPUTS REQUIRED ARE YD-HYDRIDE DENSITY, YC-HYDRIDE CAPACITY, PD-PA 

CKING DENSITY, H-HYDROGEN TO BE ABSORBED 

50 REM INPUTS REQUIRED ARE CD(1)-MAX COIL DIA, CD (2) -MINCOIL DIA, SH-MAX S 

TACK HEIGHT, OD-TUBE OD, WT-TUBE WALLTHKNS, SP-SPACE BETWN COILS, FD-AX IAL FILTER 
DIA, SL-SPACE BETWEEN LAYERS 

60 REM OUTPUTS ARE TL (3) -TOTAL COIL LENGTH PER LAYER, TT-TQTAL COIL LENGTH 

PER STACK, YW-HYDRIDE WEIGHT, SN-NUMBER OF STACKS REQD, NL-LAYERS PER STACK 
70 INPUT "INPUT DIAMETER OF LARGEST COIL " ; CD (1 > : INPUT " 

INPUT DIAMETEP OF SMA 
LLEST COIL ";CD<2> 

80 INPUT " 

INPUT MAXIMUM STACK HEIGHT, FT " ; SH: INPUT " 

INPUT SPACE BETWEEN CO I 
LS, IN "*SP 
90 INPUT " 

INPUT SPACE BETWEEN LAYERS " ; SL 
100 INPUT " 

INPUT TUBE OUTSIDE DIAMETER, IN "*0D: INPUT " 

INPUT TUBE WALL THICK 
NESS, IN " , WT 
110 INPUT " 

INPUT AXIAL FILTER OD, IN " 5 FD 
120 INPUT " 

INPUT HYDROGEN TO BE STORED, LBS "*H<3) 

130 ID - OD - 2 • WT:Y « 1 

140 P » OD ♦ SP:NC<1) • ( CD ( 1 ) / 2) / (P / 12):NC(2) - <CD<2> / 2) / (P / 12) 

150 TH < 1 ) - NC(1) • 2 • 3. 14: TH <2> - NC(2) • 2 • 3, 14 

160 A * CD < 1 > / (2 • TH ( 1 ) ) : TL < 1 ) - A * TH <1 ) s 2 / 2:TL(2) - A • TH<2) ~ 2 / 2: 

TL(3> - TL < 1 ) - TL (2) 

170 INPUT * 

INPUT HYDRIDE DENSITY, LB/FT3 " « YD: INPUT " 

INPUT HYDRIDE CAPACITY, 

7. " » YC: INPUT " 

INPUT PACKING DENSITY, 7. *' s PD 

ISO NC <3> - NC(1> - NC ( 2 > : X - NC(3>: GOSUB 780:NC(3) - X: PRINT '* 

TOTAL NUMBER 

OF COILS PER LAYER - ,, ;NC<3) 

190 PRINT D**"PR#1": PRINT I**"U"I PRINT I*| "SON" 

200 ID - OD - 2 * wi 

210 V - ((ID ~ 2 - FD 2) • .7054 * TL(3>> / 144: YC - YC / 100: PD - PD / 100 

220 YW<3> - H(3) / YC:YW(1> - V • PD * YD: H ( 1 ) - YW ( X ) * YC 

230 X - TL<3>: GOSUB 700 

240 TH3) - x:x - V: GOSUB 700 

250 V - x:x - HY: GOSUB 700 

260 HY - x:x - H: GOSUB 700 


270 

PRINT 

'• •• 






200 

PRINT 

TAB ( 

25) * 


INPUT 

DATA *♦♦***" 


290 

PRINT 

" 






300 

PRINT 

TAB ( 

5) \ " 

TUBE 

DIAMETER - 

" s OD* " IN"? TAB < 15)* "WALL THICKNESS - 

"1 WT 

IN" 







310 

PRINT 

M .1 





" ? CD 

320 

PRINT 

TAB ( 

5) ; « 

'COIL 

DIAMETER — MAXIMUM - " 5 CD < 1 ) ? " FT MINIMUM - 

(2> * 

» F T " 








330 PRINT " " 

340 PRINT TAB < 5) I "SPACE BETWEEN COILS « "JSPI" IN"; TAB ( 12) ; "BETWEEN LAYERS 
- "tSLI" IN" 

350 PRINT " " 

360 PRINT TAB < 5) * "HYDRIDE CAPACITY - " * YC * 1001" 7." i TAB < 15) 5 "PACKING DENSI 
TY - " ; PD * 100* " 7." 

370 PRINT " " 

300 PRINT TAB ( 5)*"MAXIMUM STACK HEIGHT - "?SH*" FT" I TAB < mi" AX IAL FILTER 0 
D " * FD* " IN" 


TABLE 4: 


Computer listing of program to calculate H2 parameters of a spirally 
wrapped layered stack. (Continued on page 6.) 



390 PRINT " " : PRINT " " 

400 PRINT TAB ( 25) ; "*♦***♦ LAVER TOTALS ♦***♦*" 

410 PRINT M PRINT •• •• 

420 PRINT TAB ( 5) ; "NUMBER OF COILS * ";NC(3)5 TAB ( 19) I "COIL LENGTH - ";TL<3)i 
" FT" 

430 X - YW<1>: GOSUB 780 
440 VW(1> - X 
430 PRINT " " 

460 PRINT TAB ( 3): "HYDRIDE WEIGHT * "YW<1)|" LBS"; TAB < 11); "HYDRIDE VOLUME - 
" ; V * PD f " FT3” 

470 NL * INT (SH * 12 / <0D + 3L)> 

480 PRINT *• PRINT " PRINT " " 

490 PRINT TAB ( 25) 5 "*♦***♦ SYSTEM TOTALS ******" 

300 PRINT " "5 PRINT " " 

310 SN - (YW<3) / YW ( 1 ) ) / NL 

520 X - SN: GOSUB 780 
530 SN < 1 ) - X 

340 PRINT TAB ( 5) 5 "COIL LAYERS PER STACK - "|NLJ TAB ( 12) I "NUMBER OF STACKS RE 
QUIRED - ; SN ( 1 ) 

530 PRINT " " 

360 TT - TLC3) • NL • SN: X - TT: GOSUB 780 
370 TT - X 

380 TV * V • NL • SN:X - TV: GOSUB 700 

390 TV - X 5 TV ( Y ) - TV • PD: X - TV(Y): GOSUB 780 

600 TV ( Y ) - X: TV(K) * TV(Y) / 33-3:X - TV<K): GOSUB 7B0:TV(K) - X 

610 PRINT TAB ( 5> ; "TUBING LENGTH « INT (TT) 5 " FT"; TAB ( 14) INTERNAL VOLUM 

E - " ! TV| " FT3" 

620 PRINT " " 

630 PRINT " " 

640 PRINT TAB ( 3) i "HYDRIDE VOLUME * ";TV(Y)5" FT3 "?TV<K>5" M3" 

650 TK - YW(3> / 2. 2: X - TK: GOSUB 780: TK * X 
660 PRINT " " 

670 PRINT TAB ( 3) 5 "HYDRIDE WEIGHT - "5YW(3)|" LBS "5TKJ" KGS" 

680 HW » YW(1) * NL • SN • YC:X * HW: GOSUB 780 
690 HW - X 

700 HK - HW / 2.2: X - HK: GOSUB 780:HK - X 

710 CF - HW / . 0055 : X - CF: GOSUB 780 

720 CF - XICM - CF / 33. 3: X - CM: GOSUB 780: CM - X 

730 PRINT " " 

740 PRINT TAB ( 3) I "HYDROGEN STORED - ";HWf" LBS " $ HK ; " KGS" 

730 PRINT " " 

760 PRINT TAB ( 3) ; "HYDROGEN STORED • "jCFi" SCF ";CM;" SCM" 

770 GOTO BOO 

780 X - ( INT (X • 100 + .5)) / 100 
790 RETURN 

800 CY ■ 14 • TK:X * CY: GOSUB 780:CY * X 

810 CS « 1.4 • CY:X » CS: GOSUB. 780:CS - X: PRINT " " 

820 PRINT TAB ( 3) I "HYDRIDE COST - S"; INT (CY) J TAB < 17); "SYSTEM COST - 4"; IN 
T (CS) 

830 PB « CS / 1.8E6;X - PB: GOSUB 780: PB - X 

840 PRINT " PRINT TAB ( 5) j "PAYBACK « * 1 . 8MM/YR « "?PB; M YEARS" 

900 PRINT CHR* (12): PRINT D0| M PR#O M : END 


TABLE 4s Continued 


HYDRIDE CAPACITY, % 

1.0 

1.2 

1.4 

TUBING REQ’D, 10 5 ft 

3.40 

2.83 

2.43 

HYDRIDE REQ’D, 10 5 lb« 

2,40 

2.00 

1.77 

SYSTEM COST,' 10 6 $ 

4.33 

3.01 

3.09 

STACK HEIGHT, ft 

mi 

1 5 

20 

10 

15 

20 

10 

15 

20 

COIL DIAMETER, ft 

i i 

5 | io\iZ'l5 

B3 

BB 


Bl 

0 

10 

i 

/2|/5 

5 loi/2|/5 

I i 

S j 10 j 12 IS 

— 1 — 
5 /0|/2 

IS 

■5\I0\I2 IS 

HMj 

VESSELS REQUIRED 

zskLta 

1 1 1 ! 1 

m 



as 

s 

a 

as 

mmm 

ttnntti 

12 I 3 izi / 

II: 

20l 4 i 3 

1 I 

z 

EBBS 

io\z\iAi 

iii 

PAYBACK TIME, ymm 

2.40 

2.00 

1.72 


TABLE 5: Computer projections for typical capture systems 





































n INTRODUCTION 


Several million gallons of liquid hydrogen boiloff losses are anticipated 
annually when the Space Shuttle program reaches full stride during the late 1980's 
unless means are found to prevent or reclaim this boiloff. The subject work tested 
the concept of using metal hydride hydrogen absorption systems to rapidly capture 
low pressure (<2 psig) gaseous hydrogen for subsequent reliquefaction at a greatly 
reduced rate. 

A. Background 

In 1978, Mr. Ed Snape of Ergenics and Dr. Gerry Golub of PRC Systems 
Service began discussing the possibilities for using metal hydrides to reduce 
hydrogen loss at the Kennedy Space Center (KSC). After a period of inactivity, the 
discussion was resumed between Greg Egan of Ergenics and Dr. Golub in April 1981. 

In October 1981, Frank Lynch of Hydrogen Consultants, Inc., a subsidiary 
of Ergenics and Greg Egan traveled to KSC to meet with Jim Spears of NASA and 
others to explain hydride technology and begin a detailed discussion of possibilities 
for using the hydrogen normally lost during KSC operations. 

The first possibility, which was considered in earlier discussions between 
Snape and Golub, was to fuel the buses which carry visitors on tours of the KSC. 
This presented a number of operational problems since NASA has developed 
hydrogen handling procedures and regulations which would need many modifications 
to permit the use of hydrogen in buses. The bus conversion would also involve a 
contractor (TWA) who operates the bus fleet for NASA. 

In January of 1982, Egan and Matt Rosso of Ergenics met with Spears 
and others at KSC to present other possibilities for LH 2 boiloff utilization and 
re liquefaction. The possibilities considered were to: 

1. Use steady-state boiloff to fuel a 50 KW Motor Generator (M-G) set. 

2. Alter LH 2 delivery schedules so that transfer losses would be more 

uniform thus permitting the use of larger M-G sets producing 200-400 

KW of electrical power. 
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3. Combine M-G sets with hydride technology (storage and compression) to 

augment the performance of a reliquefaetion device. 

These suggestions were well-received by NASA personnel. After 
additional preparation, Egan, Mark Golben and Rosso of Ergenics with Lynch of HCI 
returned to present a more detailed discussion of these concepts in April 1982. 
During this meeting it became clear that one of the most attractive possibilities for 
hydride use at KSC was the capture of large boiloff losses during transfer operations 
since this would complement the reliquefaction processes currently under 
consideration by NASA. 

As the team of hydride experts became more familiar with the KSC 
Standard Operating Procedure (SOP) it became increasingly apparent that viable 
technical projects would require a sizeable effort to become familiar with the 
details of the hydrogen systems and with the people, procedures and regulations 
which control their use. 

An unsolicited proposal was submitted in June, 1982 seeking a contract 
to demonstrate the technical feasibility of the most critical part of the hydrogen 
boiloff capture/reliquefaction concept - - rapid capure of the low pressure hydrogen 
vapor. In April, 1983 the subject contract (NAS10-10625) was awarded to Ergenics 
to perform the ten month program. 

B. Hydride Applications at KSC 

Most of the hydrogen vented at LC-39 is released in rapid bursts of short 
duration while cooling transfer plumbing or relieving pressure. It is uneconomical to 
install liquefiers large enough, or to use conventional gas compression and storage 
devices, to counter these losses. Metal hydrides, however, are capable of rapidly 
absorbing hydrogen gas at low presssures and normal temperatures and storing the 
captured gas until an economically sized liquefier can return it to the storage tank 
as liquid. 


Figure 1 shows the cumulative loss estimate during one year of 
operation. At the bottom of the graph the relatively small stable boiloff is shown. 
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All of the reliquefaction schemes under consideration at NASA have the capability 
of capturing the stable boiloff so no hydrides are necessary for that purpose. 



LAUNCH AND SCRUB EVENTS DURING ONE TEAR 


Figure 1: Cumulative boiloff losses during one year. Capturable losses are shown 

as cross-hatched area. 

The region just above the stady-state boiloff line in Figure 1 represents 
the cumulative loss from truck and railcar deliveries of liquid hydrogen. Over one 
million gallons will be lost each year during these short (1.0 hour) LH 2 delivery 
periods. 


The contribution of the pre-launch, post-launch and scrub losses 
represents more than 50% of the total boiloff; but hydrogen from these sources is 
not readily recoverable since they are too intimately associated with shuttle launch 
operations. Hydrogen vapor generated during the "pre-launch” operation is the 
largest single source of vented vapor. This hydrogen is vented at the Space Shuttle 
during the cooldown and filling of its external tanks (ET's), where as much as 
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250,000 gallons are lost during the refueling operation. The "post-launch" losses 
result from blow-down of the fill line which transports liquid hydrogen from the LH 2 
storage dewar to the launch pad. The pressurized contents of this line is typically 
drained back to the dewar until pressures are equalized. The residual, low-pressure, 
hydrogen is then sent to the bum pond where it is ignited. 8000 gallons are flared 
for each launch episode which may be recoverable in the future when confidence in 
the capture system has been established. Capture of losses resulting from scrubbed 
missions, although large (100,000 gallons), are not cost effective since they are 
infrequent, unpredictable and too intimately associated with shuttle launch 
operations. 


The present feeling at the KSC is that it is preferable to vent the 
hydrogen from the above three sources rather than risk over-pressurizing the 
external liquid hydrogen tank by preventing the free flow of hydrogen boiloff. In 
this study, we, therefore, addressed only the capture of hydrogen from the two 
remaining sources; the daily, stable boiloff and the LH 2 tanker off-loading losses. 

During the transfer of hydrogen from the liquid hydrogen delivery trucks 
to the pad storage tank, approximately 10% of the fuel is lost. Five trucks arrive at 
the site, each with 13,000 gallons of LH2» The hydrogen boiloff available for 
capture was estimated by NASA personnel to be about 6,500 gallons. "First-hand" 
observation of the off-loading procedure by Ergenics personnel and a fundamental 
modeling analysis conducted during the course of this program indicate that only 
6,000 gallons are lost during this operation. 

It was assumed, as shown in Figure 2, that the largest fraction of this 
hydrogen was boiloff produced during cooling of the transfer piping. In reality, the 
results of the study portion of this program show that the greatest contribution to 
boiloff is made by flashing of the liquid as the pressure decreases from the tanker 
transfer pressure (20-25 psig) to the dewar storage pressure (2.5 psig). See Appendix 
a( 2). The time scale of Figure 2 is based on the NASA estimate of 1.5 hour total 
delivery period, with 1/2 hour of actual LH 2 transfer to the pad storage tank. Two 
"off-loadings" witnessed by Ergenics personnel required only one hour with LH 2 
being transferred to the dewar for the entire period. The boiloff rate was more 
controlled than was originally anticipated with vapor being vented to the 
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atmosphere at a fairly uniform rate. The expected initial rapid discharge of vapor, 
as the transfer line is cooled to LH 2 temperatures, was not evident; and based upon 
the modeling study is not a significant contributor to the total boiloff. The dotted 
line on Figure 2 is more representative of the actual boiloff rate. 



FIGURE fc Estimated boiloff profile during a tank fill event. 

C. Capture Demonstration 

Capture of this vapor is not readily accomplished bv normal means 
because of the extremely high mass flow-rate (>10,000 SCFM) and low presure (<2.5 
psig). Certain low pressure metal hydrides operating near ambient temperatures can 
absorb hydrogen very rapidly and well below atmospheric pressure. In the course of 
this work a proof-of-concept and a field test device were fabricated to demonstrate 
the rapid absorption of low pressure hydrogen. The body of this report contains a 
description of the device, the test results and recommendations for additional work 
culminating in the implementation of a full-scale capture/reliquefaction system at 
the KSC. 


HI. APPROACH 


The approach is described in Section 4.0 TECHNICAL REQUIREMENTS 
section of the original Statement of Work, which follows; 

4.0 TECHNICAL REQUIREMENTS 

4.1 The effort will include the following tasks: 

4.1.1 Task 1 - System Performance Requirement 

The project team will meet with KSC personnel in order to study 
and specify the following boiloff capture design criteria; 

4. 1.1.1 Pad Storage Tank Fill Operation Boiloff: 

Pressure 

Temperature 

Flow 

Time Allotment 

4. 1.1.2 Hydride Capture Unit to Liquefier Transfer: 

Pressure 

Flow 

Time Allotment 
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4.1.2 Task 2 - Hydride Alloy Selection 

Based on the results of Task 1, one or more hydride alloys will be 
specified. Preliminary estimates of thermal ballast fraction, heat 
capacity, thermal conductivity, and heats of formation will be 
made. 

4.1.3 Task 3 - Hydride Container Design 

The hydride alloy selected in Task 2 and the performance 
requirements of Task 1 will be used as inputs to the design of a 
suitable container which must provide for ample gas flow means 
for heating and/or cooling and for control of powder expansion. 

4.1.4 Task 4 - Alloy/Container Design Verification 

A small batch of one or more alloys selected in Task 2 will be 
formed and inserted in a small-scale container representative of 
the design created in Task 3. The assembled unit will be tested on 
a Sieverts apparatus to verify performance. 

4.1.5 Task 5 - Hydride Alloy Manufacture 

The selected hydride alloy, verified in Task 4, will be formed in a 
quantity sufficient for the construction of a proof-of-concept 
boiloff capture unit. 

4.1.6 Task 6 - Container Construction and System Assembly 

The basic container design verified in Task 4, will be followed in 
constructing the proof-of-concept boiloff capture unit. The 
hydride alloy will be inserted, and readied for performance testing. 
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4.1.7 Task 7 - System Test 


The completed boiloff capture unit will undergo a series of tests 
representative of the KSC pad storage tank transfer operations, 
specified in Task 1. 

4.1.8 Task 8 - Report Preparation 

A final project report will be prepared, detailing the design and 
performance of the boiloff capture unit. The report will also 
include the proposed future work statement and discussion of a 
full-scale system design for LC-39. 

4.1.9 Task 9 - Proof-of-Concept Demonstration and Project Review 

The complete boiloff capture unit will be shipped to KSC and 
demonstrated to key personnel. While the demonstration is being 
performed, it will be reviewed by the NASA project team who will 
notify interested parties at KSC who want to witness the 
demonstration. 

The task sequence will be used as a format for the technical discussion. 

A. Task 1 - System Performance Requirements 

It became apparent during discussions with NASA and EG<5cG personnel 
that very little "hard” data existed concerning hydrogen vapor mass-flow-rates. It 
was also recognized that to attempt to make definitive measurements would be a 
major undertaking. We, therefore, agreed to base the engineering design of the test 
unit on whatever subjective information could be compiled as well as projections 
generated by computer modeling from first principles. Two liquid off-loadings were 
viewed by Ergenics personnel and some subjective measurements made. The 
estimated vapor flow-rates, based on these data, agreed very well with projections 
made by independent consultants working from engineering drawings of the LC39A 
LH2 Storage Area (See Appendix A). The drawings and whatever other data were 
available were supplied by the NASA Technical Representative. The peak and 
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average hydrogen flow-rates were taken to be 10,650 and 10,300 SCFM (58.6 and 
56.6 Lb/Min) respectively. It was also determined by thermodynamic analysis that 
the temperature of the gas entering the recovery unit will be from 35 to 85° K for 
capture line lengths up to 500 feet long. The relatively small temperature increase 
is attributed to the high linear velocity of the gas and, consequently, short residence 
time. The sensible heat capacity of the cold vapor will assist in the removal of the 
exothermic heat of reaction which is generated as the hydrogen is absorbed by the 
hydride. 


In practice, the delivery trucks undergo a pressure building step in order 
to force the liquid through the transfer plumbing. The pressure within the truck's 
dewar is increased to 20-25 psig by vaporizing some of the liquid in an external heat 
exchanger. The pressure is maintained on the liquid until the truck's tank has 
discharged all of its LH 2 . The first liquid to enter the transfer line plumbing boils 
away until cryogenic temperatures are reached. 

After the transfer line has been chilled to allow liquid flow, the second 
opportunity for hydrogen capture occurs. As the pad storage tank is filled, the cold, 
gaseous hydrogen in the tank's ullage is displaced and vented. This gas is available 
for capture during the entire 1 hour fill period. After transfer, the high pressure 
gaseous hydrogen which remains in the trailers is vented before the trucks depart 
and is, therefore, available for capture during the latter stage of the delivery cycle. 

The shape of the LH 2 boiloff profile, shown in Figure 2, is quite similar 
to the shape of a hydride charging profile, as shown in Figure 3. Both curves are 
shown together in order to illustrate their similar characteristics of initially high 
flowrates followed by a period of lower flowrates. The hydride charging profile is 
taken from laboratory data for a hydride with an absorption pressure of 0.5 
atmospheres absolute at 25° C and an applied hydrogen pressure of 1 atmosphere 
absolute. The dotted line represents the more representative hydrogen boiloff rate 
as measured during this work. The actual capture mission is less demanding than 
originally anticipated, as the slope of the dotted curve indicates. 
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FIGURE 3: Comparison between the estimated boiloff loss profile and a hydride 
charging profile shows the feasibility of boiloff capture. 


The parameters used for the design of the full-scale recovery system 


follow; 


Hydrogen flow-rate (PEAK) 10,600 SCFM 
Hydrogen storage capacity 3400 lbs. (620,000 SCF) 


C’ »xl + Am / 

uyiu naici aiui c/ 

flow-rate 20° C(68° F)/400 GP M for a 1 0° C (1 8° F) Rise 

Hot water temperature/ 

flow-rate 80°C(176°F)/40 GMP for a 10°C (18°F) Rise 

Maximum vessel diameter/ height 15 ft/20 ft 

Tubing diameter/ wall thickness 1.125 in/0.050 in 

Axial filter OD 0.225 in 

hydride packing density 60% 
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B. Task 2 - Hydride Alloy Selection 


The critical criteria that impact alloy selection will be discussed for 
each alloy considered. They include: 

. Hydrogen absorption capacity 
. Chemical kinetics 

. Dynamic hysteresis 

. Isotherm slope 

. Dynamic absorption plateau pressure <0.4 atma at 25° C 

. Alloy components availability 

. Dynamic desporption plateau pressure >2 atma at 75° C 
DEFINITIONS 


Hydrogen Absorption Capacity is the amount of hydrogen that can be 
reversibly stored in a metal hydride and is usually given as weight per 
cent of the metal, (W^/WHydride) x 100. 

Chemical Kinetics is a measure of the time required for a hydrogen 
molecule to attach to the metal surface, dissociate into two H+ ions and 
be absorbed into the crystal lattice of the metal. Ideally the chemical 
kinetic rate is measured independently from heat transfer effects; but 
can have a measurable influence on the effective hydrogen absorption 
rate. Figure 4 depicts the absorption rate to 50% of full capacity for 
three different alloys. Differences greater than an order of magnitude 
are attributed to the inherent chemical kinetics for each alloy. These 
curves were generated at Ergenics under contract to Sandia National 
Laboratories in Albuquerque, New Mexico(3). Figures 5-9, also taken 
from the Sandia final report, present the measured wide range of 
absorption rates for some of the alloys considered. Figures 5 and 6 
depict the typical variation in kinetic rates as a function of temperature. 
Even small additions of aluminum seem to slow the chemical kinetic rate 
and increase its sensitivity to temperature change. This result makes it 
imperative that we know the temperature of the hydrogen vapor entering 
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the boiloff capture vessel. The effect of alloying elements, other than 
aluminum, is evident when Figures 7-9 are compared. Data are not 
available for LaNi4.gAlo.45 but reliable extrapolations can be made from 
the information presented here. 


LaNI 3 Co 2 


P .<0.35 atma 
A 


LaN, 4.9 AI 0.1 
Ss *s s P a * 1.63 atma 


P *2.31 atma 


LOG (P/P A ) 


FIGURE 4; Absorption rate to 50 % of full hydrogen capacity for three AB5 alloys. 
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FIGURE 5; Comparison based on reduced pressure of LaNig absorption kinetic data 
showing the rate and time to 0.05 g-atom H reacted at all pressures and 
temperatures tested (^0.5 H/M). Ref. 3. 



FIGURE 6; Summary of absorption kinetic data at 25° C showing the maximum 
temperature excursion as well as the reaction rate and time at 0.05 g- 
atoms hydrogen reacted (^0.5H/M). Ref. 3. 
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FIGURE 7: Summary of absorption kinetic data at 25° C showing the maximum 
temperature excursion as well as the reaction rate and time at 0.05 g- 
atoms hydrogen reacted (u'O.SH/M). Ref. 3. 
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FIGURE 8s Summary of LaNi3C02 absorption kinetic data at 25°C showing the 
maximum temperature excursion as well as the reaction rate and time 
parameters. Ref. 3. 
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FIGURE 9: Comparison based on reduced pressure of LaNi^gAlo.i absorption 
kinetic data showing the rate and time to react 0.05 g-atom H at all 
pressures and temperatures tested (>r0.5 H/M) Ref. 3. 

. Dynamic Hysteresis - is the difference between the absorption and the 
desorption plateau pressure for a particular alloy while in the transient 
mode. The plateau pressures are significantly different for different 
hydrogen mass flow-rates. The desorption and absorption plateaus are 
depressed and raised respectively when compared with its equilibrium 
isotherm, as can be seen in Figure 10. For applications, such as the one 
being considered here, where hydrogen flow-rates are extremely rapid, 
the system design must be based upon dynamic values. Because the 
magnitude of this loss is an inherent characteristic of each hydride 
formulation, it is important to consider it when evaluating alloys for any 
application. 

. Isotherm Slope - some hydrides, such as LaNi3.gFei.2 (Figure 11) and 
LaNi 4 . 7 Sno .3 (Figure 12), exhibit a sloping plateau which has a profound 
effect on a hydride's applicability. For example, the tin substituted 
hydride (Figure 12) has a plateau that slopes from 0.15 to 1.5 atma at 
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. Dynamic Absorption Plateau Pressure - is a primary aspect to be 

considered since the greater the pressure difference between the 

charging hydrogen stream and the hydride's absorption plateau, the 
greater will be the absorption rate. We, therefore, want this value to be 
as low as possible at normal temperatures (<0.4 atma at 25° C). 

. Availability of Alloy Components - while not a major concern, should be 
addressed when making a hydride selection. The use of rare, expensive 
or sensitive materials must be avoided particularly where large 
quantities of hydride are to be used. Of the five candidates considered 
only the cobalt containing alloy might fit in this category. 

. Dynamic Desporption Plateau Pressure - must be high enough to 

generate pressures above atmospheric pressure at moderate 

temperatures. Liberating the captured hydrogen with the use of low- 
grade waste heat from some other process is extremely attractive from 
the economic point-of-view. Hydrogen release can be considered an 
extreme bottoming cycle for any co or poly generation system (target Pd 
>2 atma at 80®C). 

Initially five alloy compositions were considered for this application. 
Four of the candidate alloys were lanthanum-nickel with small addition of various 

other metals (Al, Co, Fe, Sn). The fifth alloy considered contained Ca, Ni and AL 

See Figures 11-15. 

Alloys I and n, the iron and tin substitutions, LaNi3.3Fei.2 and 
LaNi 4 # 7 Sno .3 respectively, were eliminated because of their sloping plateau and 
greatly reduced capacity, See Figures 11 and 12. Even though the plateaus were 
approximately at the correct presssure, demonstrated poor kinetics was also a 
contributing factor in this decision. 

Alloy HI (CaNi^gAlo.i), although it had a flat plateau at the right 
pressure, its excessive hysteresis and low capacity caused it to be rejected., The 
undesirable hysteresis characteristic would require water temperature to be greater 
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than 80°C to raise the desorption plateau above the desired discharge pressure of 2 
atmospheres absolute, see Figure 13. 
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FIGURE 13; Dynamic pressure-composition isotherms for CaNi^gAloj. Ref. 4. 

Candidate alloy IV, (LaNi3Co2>, or a slightly modified version of it, will 
probably work very well; but because it contained cobalt which might be a problem 
to obtain in large quantities, it too was eliminated. Figure 14 depicts its acceptable 
properties and Figures 4 and 8 show its moderately poor kinetics. 
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FIGURE 14s Dynamic pressure-composition isotherms for LaNi 3 Co 2 « Ref. 4. 

Alloy V (LaNi^yAlo^-HYSTOR 207) closely approximated all the desired 
qualities except that the absorption plateau pressure was too high. By comparing 
Figures 15, 16, 17 the negative effect on the hydrogen storage capacity of aluminum 
substitutions for nickel is apparent. Greater additions of aluminum would have two 
major effects on the alloy developed; one desirable and one undesirable. The 
absorption plateau would be further depressed to the desired pressure but further 
reduction in hydrogen capacity would undoubtedly result. Because past work in this 
area had shown that the capacity loss would be small; and because the kinetics of 
this family of alloys is known to be extremely fast, we chose to melt a sixth alloy 
with an even higher aluminum content. 

Alloy VI (LaNi 4 > gAlo .4 Heat T-88855-2) was melted and annealed 
(1125°C for 17 hours) and although the exact composition was not attained, the 
plateau depression can be seen by comparing Figures 15 and 18. 
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FIGURE IS: Dynamic pressure-composition isotherms for FIGURE 16: Dynamic pressure-composition isotherms for 

LaN *4.7 A lQ.3* Ref. 4. LaNi^gAloj. Ref. 4. 
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When the natural log of pressure is plotted vs. reciprocal absolute 
temperature, for a fixed hydrogen concentration (usually mid-plateau), the result is 
linear and is referred to as a Van't Hoff plot. Figure 19 is a composite plot of a few 
of the alloys considered. LaNis is also presented for reference purposes. It can be 
seen that Alloy VI can absorb H 2 at sub-atmospheric pressures at 25° C and can 
liberate the gas at pressures >1 atmosphere at 80° C. See Table 6 for the properties 
of the candidate alloys. 
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FIGURE 19 ; Van't Hoff diagrams for candidate alloys with LaN^ for reference. 
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C. Task 3 - Hydride Container Design 


The absorption of hydrogen by metal hydride is a highly exothermic 
reaction and conversely is endothermic upon desorption. The rate at which the 
hydrogen must be absorbed or desorbed determines the method of heat management 
used. Two approaches generally used in the design of hydride vessels is an 
isothermal design, where the heat of reaction is either removed or added by means of 
a heat exchanger with a thermal transport fluid and an adiabatic design where the 
heats of reaction are stored within the bed by the incorporation of large quantities of 
high heat capacity-inert materials. Since the heat transport paths are shorter in the 
adiabatic system, that approach is used for the higher hydrogen 
absorption/desorption rates. Where hydrogen flow-rates are low enough the 
isothermal process is used since smaller and fewer vessels are required with a greatly 
reduced hydride inventory leading to a desirable effect on the economics of the 
system. 


Because of the encouraging results obtained during the performance of 
Task 1, a decision was made to use an isothermal design for the management of the 
generated heat. The rapid hydrogen flow-rate made it mandatory to keep the 
hydride "bed" thin thus reducing the length of the heat transfer path. During 
absorption, cooling water is made to flow around small diameter (.375 in. OD) tubes 
in which the hydride is contained. A flexible filter, of cylindrical cross section, is 
placed in the tube with the hydride occupying the annular space between the filter 
and the cylinder wall. Since the filter runs the full longitudinal length of the tube it 
provides a low pressure-drop path for the low pressure hydrogen vapor to enter the 
system. 


The gas is very rapidly absorbed by the metal to which it is exposed and 
the heat of reaction dissipates through the outer wall into the cooling water. For 
details of construction see Figure 20. As shown on the same figure, four of the 25 
ft. long 0.375 in. OD tubes are placed into a 1.125 in. OD water jacket, and spirally 
wound to a 27 in. diameter flat coil. Three coil assemblies were fashioned in this 
way, placed inside an enclosure and the water jackets connected in series. The 
hydrogen lines, each equipped with a filter, pressure relief valve and a manual shut- 
off valve were connected in parallel to a central manifold through which the 
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hydrogen enters and exits the system. The enclosure also housed two electrically 
operated solenoid valves to control the hot (for desorption) and cold (for absorption) 
water flow through the jacket. Two motor-operated rotary-valves provide a low- 
pressure path for the hydrogen to flow from the source to the ’’bed”. A differential 
pressure transducer, located between the two rotary valves, enabled monitoring of 
the pressure in the coil, at the source or both. The transducer range, -10 psi to +10 
psi, accurately transmitted the low pressures and small changes witnessed during the 
operation of the device. A thermocouple in the hydride, another in the water jacket 
and a third in the air enabled recording of hydride, water and ambient temperatures 
during a test run. A pressure transducer (0-600 psig) was installed on the high- 
pressure side of a regulator to track pressure changes in the hydrogen cylinder 
source as the hydrogen was absorbed by the hydride. 


i 
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All temperatures and pressures as well as all electrical valve switches 
were monitored and operated remotely, from 50 feet away, for safety 
considerations. 

D. Task 4 - Alloy/Container Design Verification Test 

The design described in section C, i.e. 0.375 in. diameter copper tube 
with axial flexible filter and hydride filling the annular space between the two, had 
to be verified by flow tests. A single coil 10 ft. long was fabricated using 380 grams 
(38g/ft) of the specially melted alloy (LaNLj.gAlo^ heat T-88860-2) after grinding 
to -35 mesh. The hydride was vibrated into the tube/filter assembly’s hydride space 
where the resulting void fraction is typically 40%. The tube was equipped with a 
pressure relief valve and a manual shut-off valve. The coil was not fitted with a 
water jacket; but was immersed in an agitated temperature-controlled bath instead. 
This provided better temperature control and nearly isothermal test conditions. 

The verification test program emphasized the absorption half of the 
cycle since the rapid absorption of the low pressure hydrogen vapor was recognized 
as the constraining factor in this application. Hydrogen gas at various pressures 
slightly above and below atmospheric pressure was used to charge this coil while it 
was held at 25° C by the water bath. The time to reach 90% of full charge was 
recorded and the results are presented in Figure 21 and Table 7. At the system 
design pressure point (2 psig) the coil exhibited the capability to absorb 44 liters of 
hydrogen, 90% of the hydride's full capacity (1.16 wt%) within 2 minutes. This rate 
is 30 times faster than required for the subject application. The results of the tests 
verified the present heat exchanger design and suggested newer and less expensive 
approaches to hydride encapsulation. 
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Hydrogen 

Pressure 

psia 

Time 

To 

90%, Min 

Abs. Rate 
Max, SLPM 

Abs. Rate 
Avr, SLPM 

7 

>7.35 

15.0 

6.0 

20 

4.35 

18.4 

10.1 

15 

2.58 

32.2 

17.1 

20 

1.25 

56.0 

35.2 

25 

1.15 

80.0 

38.3 

*17 

2.0 

40.0 

22.0 


* Interpolated values for dewar design pressure. 

TABLE 7. Hydride Test Coil absorption rates at various charging pressures 
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FIGURE 21; 25° C absorption rate performance of NASA test coil as a function of 
charging pressure. 


- 35 - 


E. Task 5 - Hydride Alloy Manufacture 


Based upon the system design constraints and known hydride 
characteristics, a custom alloy of Lanthanum Nickel Aluminum was identified as the 
"hydride-of-choice". As described in Section B, Task 2, the hydride contained a 
higher aluminum content than Ergenics standard HYSTOR-207 (LaNi4.7Alo.3) alloy. 
A 30 lb melt was made in a vacuum induction furnace at the International Nickel 
Research and Development Center (IRDC) in Sterling Forest, NY. The melt was 
supervised by Dr. G. D. Sandrock of the Ergenics staff and the Chemical Analysis is 
presented in Table 8. The billet was normalized by heat treating at 1125°C for 17 
hours and crushed in a Gyro Mill grinder to -35 mesh. A small representative sample 
(8.0 grams) was enclosed in a hydriding reactor and both absorption and desorption 
equilibrium isotherms were generated at three different temperatures. Figures 22- 
23 depict the measured Pressure-Composition-Temperature (PCT) characteristics of 
heat T-88860-2 melted for the application. 


Analysis Heat Constituents, Weight Percent 

Number Date Number LA Ni A1 O2. _N2_ _C 


35631 8/12 T-88855 

Target 33 65 2 .05 .01 .02 

Actual 33.2 65.1 1.7 .02 .01 

Atomic Formula - Lai.o2Ni4.73Alo. 2 


35713 8/16 T-88860 

Target 33.1 64.3 2.6 .02 .01 .02 

Actual 32.9 65.0 2.1 .01 .002 

Atomic Formula - La1.00Ni4.67Al0.33 


I 


i 

Table 8; Chemical Analysis Report 
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FIGURE 23: Equilibrium (Static) and Transient (Dynamic) 25° C Absorption and 
Desorption Van't Hoff curves for LaNi 4 .gAlo .4 T-88860-2 
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An earlier melt, heat T-88855-2, did not meet composition or 
performance specifications. See Table 8 and Figure 18. Figure 22 shows the 
Pressure-Temperature relationship with the hydrogen concentration held constant at 
the isotherm plateau mid-point (H/M = 0.5). The dynamic absorption and desorption 
lines indicate that 25° C, 1.0 atma absorption and 75° C, 3.0 atma desorption is 
possible with this hydride. 

F. Task 6 - Container Construction and System Assembly 

The system design, as described in section C, is based on coiled hydride 
tubes in a water jacket, See Figure 20. Twelve 25 ft. long, 3/8 inch diameter copper 
tubes were each fitted with a 0.225 inch diameter hollow flexible filter. The tubes 
with full length filters, were attached to a vibrator and positioned to stand on end. 
The annular space between the filter and the inner tube wall was filled with -35 
mesh hydride powder to a packing density of 60%. Four tubes were inserted into 
each of three 1-1/8 inch diameter copper water jackets and soldered in place. The 
water jackets, with internal hydride tubes, were coiled into a flat spiral layer of 
approximately three wraps. 

The hydride tube ends were fitted with individual 2pm filters and 
connected in parallel to a central manifold through which the hydrogen both enters 
and exits the hydride beds. Each coil layer was equipped with a 150 psi pressure 
relief valve and a manual shut-off valve. The three iayered coils were stacked and 
bound together with their hydrogen manifolds connected in parallel, thus insuring 
free access by the entering gas to all the hydride beds simultaneously. This 
configuration provides the lowest possible pressure drop loss. The coil stack was 
placed into a steel-box-housing (30"x36"xl6 n ) and the water jackets connected in 
series. Series connection of the water lines results in greater sensible heat recovery 
and higher efficiencies. This method of cooling and heating was chosen even though 
more rapid absorption can be obtained by connecting the water jackets in parallel. 
Two thermocouples were placed in the coils; one in the hydride bed and one in the 
water jacket. The thermocouples would be used during the test program to monitor 
temperature transients as hydrogen is admitted or discharged from the coil and as 
water temperatures are changed. A third thermocouple was placed in the housing to 
monitor ambient temperature. Two motorized rotary valves were used to direct the 
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hydrogen flow from the low pressure source to the coil during absorption and from 
the coil to the low pressure gas accumulator during the discharge phase of the test 
(V— 1 and V-2 in Figure 24). The hydrogen line was also fitted with an oxygen purifier 
and a molecular sieve dryer to protect the hydride surface from poisoning by the 
inadvertant admission of air during the performance of the test. 

A low pressure differential pressure transducer (-10 psi to +10 psi) is 
located at a "tee" junction between V-l and V-2, thus enabling very accurate 
pressure measurements of the hydride coil, the hydrogen source and during 
absorption, both. A pressure transducer was also installed on the high pressure side 
of the regulated hydrogen supply to monitor changes in the cylinder pressure as the 
hydrogen is absorbed by the hydride. From these measurements of gas volume 
absorbed as a function of time, absorption rates can be calculated. Two solenoid 
valves (SV-1 and SV-2 in Figure 24) control the hot and cold water to the coil water 
jackets. External to the housing, manual shut-off valves are provided for evacuating 
and purging the system during set-up as well as for choosing the preferred hydrogen 
supply for each test. The motorized and solenoid valves are operated remotely (50 
feet) by toggle switches mounted on the face of the control panel. Temperatures 
and pressures are monitored from digital readout instruments which are also 
mounted on the panel face. See Figures 25 and 26 for photographs of capture device 
and remote control panel. 
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FIGORB 25: Hydrogen boiloff capture system. 



FIGURE 26; Remote control panel for KSC demonstration. 
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G. Task 7 - System Test 


The test of the demonstration device was conducted in two parts. The 
first, Task 7, involved testing under controlled conditions at the Ergenics, Inc. 
laboratory in Wyckoff, NJ and the second, Task 9, was conducted "in-the-field", 
under less controlled conditions, at the Kennedy Space Center. 


To simulate the low-pressure high- flow- rate hydrogen boiloff stream it 
was decided to use large balloons as the hydrogen source. This seemed, at first, to 
be a very simple and more-than-adequate solution since a balloon's internal pressure 
is extremely low (measured to be 0.01 to 0.05 psig) and essentially constant for the 
entire deflation period. Hydrogen absorption rate by the hydride is graphic and can 
be recorded as the balloon's diameter, hence its volume, is measured periodically 
during the run. A problem arose as this scheme was attempted with the small test 
coil as a result of a) residual air in the baloon due to insufficient purging, b) oxygen 
and/or water vapor permeation through the skin of the balloon or c) leached from 
the balloon material by the hydrogen. Whatever the reason, the effect was to poison 
the hydride's surface and reduce its storage capacity and chemical kinetics. The 
hydride was readily restored to its full capacity by evacuation to 1 mm Hg at 80°C 
and exposure to UHP hydrogen. Three different balloon materials were tried with 
various levels of success. See Figure 27. In the as-received" condition none of the 
materials were acceptable. Washing with freon and acetone seemed to improve 
their p erf ormance; but the best results were obtained by coating the inside of 


natural-rubber 
procedure an 
demonstration 


advertising balloons with silicone vacuum grease. Using 
eight foot diameter natural-rubber balloon was prepared for 
at the KSC with the large coil system. 


this 

the 
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FIGURE 27 Effect of hydrogen stored in various elastomers on cyclic absorption 
capacity of hydride. 

The laboratory absorption tests were conducted using a compressed gas 
cylinder as the hydrogen source and low-pressure regulation to reduce the pressure 
to atmospheric. The cooling water temperature was unregulated at an average 
temperature of 12°C and was circulated through the water jacket at a flow-rate of 
6 GPM. With valv e V-l dosed and V-2 open, the coil (hydride plateau) was measured 
at the pressure transducer as 3 psia (0.2 atma). The hydrogen absorption rate was 
measured by recording the cylinder pressure as a function of time. Ideally, as 
shortly as possible after opening V-l the system pressure should have risen to near 
atmospheric pressure, which would most closely simulate the anticipated conditions 
during boiloff recovery at the KSC. This did not happen since the hydride's ability 
to absorb hydrogen, even at the sub-atmospheric charging pressures (10 psia, 0.7 
atma) was much greater than the regulator's ability to dispense hydrogen. The 
result of this test is depicted as the solid dots (Small C v Regulator on Figure 28. A 
second run with the same test conditions, but using a regulator with a larger C v > 
resulted in the curve described by the open dots (Large C v Regulator) also in Figure 
28. Even in this run, although much improvement is evident, the coil was still 
"starved" by the regulator's inability to output hydrogen at a rate rapid enough to 
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satisfy the coil. Desorption of the coil after each run under controlled conditions 
confirmed that the maximum storage capacity, 68 SCF, of hydrogen had been 
absorbed in less than 10 minutes and that 90%, 60 SCF, had been absorbed in less 
than 3 minutes. A third test was run after returning from the KSC demonstration to 
confirm that the storage capacity and kinetics are unaffected after poisoning and 
subsequent reclamation. Using a bank of gas cylinders as a low-pressure (14.7 psia, 1 
atma) surge and "feeding" the surge through a very large C v regulator it was 
possible to run the test under conditions that more closely simulated the boiloff 
capture conditions during off-loading. The results are presented in Figure 28 as the 
open square points. The rapid initial rise, 50% absorbed in less than 1.5 minutes, is 
somewhat balanced by a slowing down above 80%. We attribute this to heat transfer 
limitations; but it should be recognized that the hydride need only reach 90% in less 
than 60 minutes to capture the off-loading boiloff. 


NASA LARGE COIL 
ABSORPTION 



FIGURE 28: NASA hydrogen capture demonstration system performance. 
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H. Task 8 - Report Preparation 


The preparation of the final report was an on-going effort with sections, 
figures and tables being prepared as required for monthly Progress Letter reporting. 
At the conclusion of the test program approximately 15% of the final report had 
been completed. 

I. Task 9 - Proof of Concept Demonstration and Project Review 

The device was operated at the KSC for evaluation by the Project 
Manager and other interested NASA and EG«5cG personnel. 

Because of the uncertainties that still existed concerning the balloon 
poisoning effect, the field test procedure was planned to first demonstrate the 
absorption of pure hydrogen from a compressed gas cylinder, the second was to 
inflate the balloon with the hydrogen liberated by the hydride as hot water was 
circulated through the water jacket, and finally an attempt would be made to re- 
absorb the hydrogen from the balloon into the hydride. 

The two regulators used at Ergenics for the first two absorption tests 
were piped in parallel so that the contribution of each might be enough to 
adequately "feed" the coil. Because of a water main break, cooling water had to be 
pumped-from a water tanker. The cooling water was warmer (19-20°C) than 
optimum and once again the hydrogen flow-rate through the regulators was 
insufficient to supply the hydrogen at the hydride's absorption rate. The system 
pressure stayed below 0.5 atma for 2.5 minutes and rose to 0.8 atma when the 
regulator setting was increased to >40 psig. The total time to 90% capacity when 
the system pressure finally reached 1.0 atma was 8.0 minutes; much longer than the 
laboratory results but still very acceptable. 

The liberation of the hydrogen from the hydride bed into the balloon was 
attempted. With V-l closed and V-2 opened SV-2 was closed and SV-1 opened. This 
permitted hot water (73° C) to flow through the water jackets, raising the hydride 
temperature and thus increasing the system pressure. In about 1.0 minute, when the 
system pressure was above 2.0 atma, V-2 was opened to permit the hydrogen to 
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inflate the balloon. The balloon inflated rapidly, accepting approximately 70% of the 
hydrogen stored in the hydride within 5 minutes. Due to windy conditions, it was 
decided to abort and begin the third test quickly. The positions of SV-1 and SV-2 
were reversed causing cold water to flow through the water jackets. With both V-l 
and V-2 open, the balloon's diameter began to decrease as the hydrogen was 
reabsorbed into the bed. Due to a number of problems only half the balloon's volume 
was absorbed into the bed. The hydride had been damaged as a result of poisoning 
by oxygen, water vapor or some other substance emitted by the balloon's material. 
The apparatus was returned to the Ergenics Laboratory where the hydride was 
regenerated and additional tests run. The tests confirmed that hydrides can absorb 
low-pressure hydrogen at a rate that is adequate to capture LH 2 boiloff during the 
off-loading operation and liberate the hydrogen at pressures greater than 
atmospheric using hot water at reasonable temperatures. 

A presentation, including slides of the figures and tables in this report, 
was made describing the fabrication and test program. Copies of all presentation 
materials and computer programs used during the performance of this program will 
be delivered to NASA/KSC with this report. 

IV POLL SCALE SYSTEM 


During the off-loading operation, hydrogen vapor is liberated at an 
average rate of 10,300 SCFM. A liquefier capable of producing 40 tons/day of liquid 
would be required to handle this off-gas in "real-time". Obviously, this is not a cost- 
effective approach. Even though a full-scale metal hydride boiloff capture system 
will require miles of tubing and many pounds of hydride, this work indicates that it 
is a feasible and economically justifiable approach to reducing launch costs at the 
KSC. See Tables 1-3, 5. 


The system must be capable of capturing and storage, albeit short-term, 
3,400 lbs. (620,000 SCF) of hydrogen for subsequent liberation to a small, 
approximately 2 tons/day, liquefier. Some combination of a smaller hydride capture 
system operating in parallel with the 12 ton/day liquefier presently being considered 
by NASA might be an even more cost-effective approach. 
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Generally, the system outline in Table 2 is representative of what a 
reasonable system would look like. For all computer runs we assumed 1-1/8 inch OD 
copper tubing with 0.05 inch wall. The space between each spirally wrapped coil was 
held at 0.5 inches as was the vertical space between the stacked coil layers. The 
axial filter was held to be 0.225 inches with a hydride packing density of 60% in the 
annular space between the filter and the copper tube wall. The hydride's bulk density 
was assumed to be 420 lb/ ft ^ which is a good approximation for the alloy to be used. 
The minimum diameter for the smallest coil wrap was held constant at 2 feet. Only 
the coil layer's large diameter, the allowable stack height and the hydride's storage 
capacity were allowed to vary. The computer program allows changing any of the 
physical dimensions or constraints, and a program disc is included with this report 
which allows NASA personnel to investigate the system's sensitivity to parameter 
changes. 


A typical system of this design would require 290,000 lbs. of hydride 
housed in 200,000 feet of tubing. With coil diameters and stack heights limited to 
10 feet, we need 5 stacks of coils, each in its own cooling water vessel. The rough 
approximation of the hydride and the system cost was based solely upon the assumed 
storage capacity (e.g. 1%, 1.2% or 1.4%) and the total hydrogen stored, 3,400 lbs. 
The system cost was assumed to be 1.4 times the hydride cost ($ 20/kg, $9/lb). For 
the computer run presented in Table 2, the hydride's storage capacity was assumed 
to be 1.2%; and the resulting costs of $2.6M for hydride, $3.6M for the total system 
and a payback time of 2 years are projected. The payback time is based on 20 
launches/year and is computed as follows. It does not include the liquefier initial 
cost, operating costs or interest on capital outlay. 

20 launches x 9 off-loadings/launchx6000 gal/ off- loading = 1.08 x 10^ gal/year 

400 gal/ day normal boiloffx360 days = 144,000 gal/year 

Total boiloff reclaimed 1.22x106 gal/year or $1.84M/year at $1. 50/gal 

The capture system can be integrated into the existing facility at LC-39 
as pictured in Figure 29. The reliquefier, being small, can be located at the storage 
dewar top where its discharge can be added directly to the existing liquid. Multiple 
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capture vessel modules with appropriate valving will permit discharging of some 
while others may still be absorbing boiloff. Additional modules may be added as 
needed. 
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FIGURE 29: Artists concept of hydrogen boiloff recovery system at LC-39. Ref. 5. 
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Figure 30 shows schematically how the capture vessel, pictured in the 
upper right-hand corner, can be incorporated into the present system with a 
minimum of interruption and no change in present operating procedure. During 
installation the existing 10 inch - 150 lb. manual valve in the vent line would be 
temporarily closed, and the flanges broken, see Figures 31 and 32. The spool piece 
between the manual valve and the 10 inch - 150 lb. 2.5 psi check valve will be 
replaced by a flanged tee with the branch leading to another 10 inch - 150 lb. 
manual valve. At this point the existing manual valve can be reopened since all 
additional construction can take place downstream of the new 10 inch manual valve. 
The location of the check valve insures that the dewar cannot be subjected to 
excessive back pressure. Should the capture vessels not be able to absorb the boiloff 
vapor for some reason, the back pressure will unseat the check valve, venting the 
vapor as is presently done. 


FILL & DRAIN LINE 



FIGURE 30 : Schematic of hydrogen vapor recovery system addition to existing LH2 
fill and drain lines. Ref. 1. 
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FIGURE 31: LH 2 dewar piping layout. Ref. 
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PIGUltB 32; Detail of hydrogen recovery system connection point in existing piping. Ref. 



V. FINAL COMMENTS AND RECOMMENDATIONS 


The program accomplished its stated mission; to demonstrate the 
technical feasibility of capturing low pressure hydrogen vapor as a reversible metal 
hydride at very rapid rates for subsequent reliquefaction. Based upon the results 
presented in this report, Ergenics believes the hydrogen capture system can be 
scaled directly for full scale operation. The extremely rapid absorption capability 
of metal hydrides suggests the alternative possibility of using a much smaller system 
operating on a short (3 to 5 minute) absorption/desorption cycle and acting as a 
staged compressor; the boiloff being stored as a compressed gas. Potential savings 
may be realized since at pressures as low as 1000 psi only fifteen tube trailers would 
be required to store all the gas for the slower reliquefaction process. Future work 
by Ergenics teamed with a suitable engineering construction firm should include: 

1. An engineering analysis and cost study be done for each of the 
capture/reliquefaction scenarios identified. 

2. Fabrication of a single full-scale component be completed for testing at 
the KSC on an actual LH 2 boiloff stream. 

3. Based upon the single component test results, design improvements 
should be added and a complete system fabricated at one launch 
complex. 

4. Similar systems to be constructed at other launch sites. 
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HYDROGEN RECOVERY, SIZING OF QUANTITY OF HYDROGEN VAPORIZED IN FILLING 


The boiloff or vaporization of LH 2 as it is ’’off-loaded" from the tankers 
to the main storage Dewar at KSC is a product of the following four contributions. 

1. Cooldown of the vacuum jacketed lines. 

2. Boiloff from flashing of a high pressure liquid into a low presssure 
volume. 

3. Boiloff from cooling down of dewar as liquid fills the dewar. 

4. Ullage displacement. 

Several steps are necessary in order to find the boiloff from each of the 
above contributions. 

1. Find the equivalent length of piping of the fill lines to the dewar. 

2. Determine the LH 2 flow rate from the known transfer line pressure. 

3. Calculate the heat capacity, cool down time and boiloff based on the 
fill-line length determined in Step 1. 

4. Determine the gas displaced from the dewar's ullage space as the liquid 
level rises in the dewar* 

5. Estimate the liquid evaporated during the cool-down of the dewar walls. 
The temperature differential between the top and bottom of the dewar 
was assumed to be 5° K. 

6. Calculate flashing into the dewar using the First Law of 
Thermodynamics and known states to determine the final quality of the 
liquid in the dewar. 
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The contributions are added and the total flow-rate as a function of 
liquid initially in the dewar may be determined. 

The actual length of the vacuum fill line is approximately 178 feet. Also 
in the line are 1 flexible hose, 1 elbow, 2 relief valves, and 3 globe valves. These 
fittings are equal to about 66 equivalent elbows. The equivalent length of piping is, 

therefore, approximately 700 feet. (1). Using the equivalent length of piping and a 
known pressure drop between the line pressure and the dewar, the approximate mass 
flow rate of liquid can be determined. 

The pressure drop (lb/ft2) in terms of viscosity, Reynold's number and 
mean fluid density is: 

P = fLu2 NRe 2/e p g cD 3 ( 2 ) Eql 


where: 

f is the coefficient of friction, dimensionless 
L is the equivalent length, ft. 

D is the inside diameter, 4 inch scd 5, ft. 

U is the viscosity, lbm/ft-sec. 

Nr 6 is the Reynold's number, dimensionless 
P is the mean density of the liquid, lbm/ft. 3 
g c is the gravity conversion constant, ft/sec2 

The same equation can be reduced to the following when known terms are 
substituted: 

P = 2.2xl0- 8 N Re 2/p Eq 2 

here the viscosity is 9.54x1 0 - 6ibm/ft-sec. at 20°K. 

Table I is a compilation of the results of the pressure drop calculations. 
The results were obtained by assuming a pressure drop, then guessing at a Reynold's 
number and getting the corresponding coefficient of friction. Using these numbers 


» 
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in the pressure drop equation, the calculated results were checked to see if both 
sides of the equation were equal. If not, then a new Reynold's number was guessed 
at and the results checked. The iteration was continued until the correct Reynold's 
number was obtained for the pressure drop assumed. The mass flow rate of liquid 
hydrogen could then be determined from the following equation: 

m = NrqvlA/D Eq 3 


where: 

m is the mass flow rate, lbm/sec. 

Nr 6 is the Reynold's number, dimensionless 
p is the viscosity, lbm/ft-sec. 

A is the internal cross-sectional area of the pipe, ft. 2 
D is the internal diameter, ft. 


Which is actually a rearrangement of the equation for computing the Reynold's 
number 


Nr 6 =rnD/»uA and as A = ttD2/4 
Nr 6 = 4m/iryD 


The results of the pressure drop and mass flow rate calculation are given 
in Table I. The equivalent volume flow rate of the liquid if all liquid were converted 
to gas at standard conditions is given in column 4. Columns 5, 6 and 7 present the 
mean density of the liquid, the velocity of the liquid in the pipe and corresponding 
Mach number, respectively. 


p 

PSIG 

N Re 

m 

lbm/sec 

Q 

SCFM 

m 

lbm/cu.ft. 

V 

ft/ sec 

M 

30 

6.8E6 

18.36 

2.11E5 

4.17 

86.6 

.077 

25 

6.2E6 

16.74 

1.93E5 

4.21 

78.2 

.070 

20 

5.5E6 

14.85 

1.71E5 

4.23 

69.0 

.062 

1 5 

4.7E6 

12.69 

1.46E5 

4.25 

58.7 

.052 

10 

3.8E6 

10.26 

1.18E5 

4.27 

47.3 

.042 

7 

3.2E6 

8.64 

9.95E4 

4.27 

39.8 

.035 

5 

2.6E6 

7.02 

8.09E4 

4.29 

32.2 

.028 


TABLE I: Mass flow rate of liquid in the fill line as a function of differential 

pressure between the fill line and dewar. 
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The next thing to calculate is the cooldown of the vacuum jacketed 
piping. There are two contributions to the liquid boiloff as it first flows down the 
piping. One term is the steady state heat transfer through the insulation. Since the 
temperature varies from approximately ambient down to the liquid temperature, 
one-half of the steady state heat gain is used. The allowable heat gain by NASA 
specifications is 1000 BTU/hr for the dewar and 1800 BTU/hr for fill lines. (2) The 
external heat transfer gain is then about 1400 BTU/hr. Assuming that the cooldown 
time is about 5 minutes, the amount of heat transferred in that time is 1.2x10^ 
Joules. This term will be shown to be negligible as compared to the sensible heat of 
the lines, which is the other contributing term to liquid boiloff during the cooldown 
of the transfer lines. The vacuum jacketed inner line is constructed of Schedule 5, 4 
inch Invar pipe. The heat capacity for nickel as a function of temperature was 
obtained from the Handbook of Chemistry and Physics, 41st Ed. page 2273 (5). The 
heat capacity for iron is roughly the same as that for nickel. The data were fitted 
to a curve and over the range of 20-300K, the equation is: 

Cp(T) = .0145+8.557x1 0~4 t - 1.979x1 0-6 t2+1.536x10"9t 3 cal/g-K 

The average sensible heat capacity over the temperature range of 20-300 K is 

0.070 cal/g-K. 

The length of tubing is about 178 feet and the mass of metal per foot is 
3.915 lbm/ft for Schedule 5, 4 inch pipe. Allowing 25% extra for supports, valve 
bodies and getter material, the mass used to determine the heat capacity was 870 
lbm. The amount of energy needed to cool down this mass is 3.27xl0 7 Joules. This 
amount is more than a factor of one hundred greater than the heat leak from 
ambient. Using the calculated heat capacity and the external heat gain, it would 
take about 2.8 days for the transfer lines to warm up to ambient temperature. 

A calculation must be done to find the mass of liquid required to cool 
down the fill line and also find the time required to cool the line down. To find the 
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time to cool down the line, it is assumed that the boiloff is so fast that choked flow 
is reached in the line, which is normally the case. The calculation is based on a 4 
inch line with a discharge coefficient of 0.60 and a cross-sectional area of 0.020 sq. 
ft. at the exit of the fill line. The equation for the velocity of the choked flow is: 

V = C d ( Ygc RTg)l/2 


where: 

V is the velocity, m/sec. 

C d is the discharge coefficient, dimensionless 
y is the ratio of Cp/C v , dimensionless 
g c is the gravity conversion coefficient, m/sec. 2 
R is the Gas Constant for hydrogen, m-Kg/Kg-°K 
To- is the average gas temperature, °K 

o 

At 1.17 atmospheres pressure and 160 K the density of hydrogen is 16.77xl0“ 3 g/cc. 
The mass flow rate of hydrogen gas is therefore .197 lbm/sec or 5.83x103 gm/min at 
choked flow conditions. This is equivalent to 2.27x1 0 3 SCFM. The time that this 
flow is taking place is only for the amount of time that it takes to cool down the fill 
lines. Time to cool down the line is given by the following equation: 

tss = (2m(Ti-T ss )Cp-V(p ss (u S s“Hf) = Pi(ui“Hf)/(mg(Hg-Hf)-0.5Qss) 
where: 

t^ is the time to reach steady state, sec. 
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Zm(Tj-Tgg)Cp is the energy to cool down the fill line due to heat capacity, Joules 
V is the volume of the line, cc 

Pss and are the steady state and initial densities in the line, g/cc 

u ss and uj are the steady state and initial internal energies of the hydrogen, J/g 

Hf is the enthalpy of the hydrogen liquid at steady state, J/g 

m<r is the mass flow rate of hydrogen gas calculated at choked flow, g/sec 

O 

Hg is the enthalpy calculated for the hydrogen gas at the average temperature in 
the line, J/g 

Qss is the steady state heat gain, cal/ sec. 

It was assumed that the quality did not change very much for the liquid 
and was considered to be zero. Doing this, the states necessary to complete the 
above equation are: 

Pss = .0663 g/cc at 2.5 atm. sat. liquid 

Hf = =219.6 J/g at 2.5 atm. sat. liquid 

Uss = -223.4 J/g at 2.5 atm. sat. liquid 

Pi = 8.1 8x1 0 -5 g/cc at 1 atm. and 300 K 
uj = 2.989x1 0 3 J/g at 2.5 atm. and 160 K 

Hg = 2. 01 2x1 0 3 J/g at 2.5 atm. and 160 K 

Substituting these into the equation for time to reach steady state, the result is 
1.47 minutes. 
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One more equation is needed and that is the mass of liquid hydrogen 
needed to cool down the fill piping. It is as follows: 


Mf/M w = C w (Ti-Tss)/(Hg-Hf)+V (p ss -Pi)/M w +Q ss t S s/ (2 M w (Hg-Hf) 


where: 

Mf/M w is the lbm of liquid hydrogen needed per lbm of metal 

The rest of the terms are as above. Substituting the necessary values the mass of 
hydrogen needed per lbm of metal is .209 lbm H 2 / lbm metal. 

The pressurization of the dewar can be calculated from the saturation 
rule given on page 475 of Barron. Thus rule, especially for hydrogen is in error on 
the liberal side, but for lack of information, this rule will serve well enough since 
the time to pressurize is on the order of only a few minutes. The saturation rule is 
the difference in the volume times density at the initial and final states. 


Mg-Vg2 Pg2 " v gl Pgl 

where mg is the change in mass of the gas from state 1 to 2 

Vg 2 and Vgi are the volumes the gas occupies at states 1 and 2 respectively 

Pg2 and pgi are the densities of the gas at states 1 and 2 respectively 

In terms of the number of gallons of the system, the volumes can be written as 
follows: 

Vg2 = 1.13xl(P _ 0.134V 

Vgl = 1.13x105 - (V + 55,000) 0.134 

where V is the initial volume of liquid in the dewar 

The saturation rule can then be written as 

m g = (1.13x105 - 0.134V)p g2 -(1.13xlo5-0.134(V+55,000))p gl 
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The time to pressurize can be determined by the maximum flow rate of 
liquid into the dewar as a function of the pressure drop. This is where much of the 
error comes in. The liquid must be able to flash entirely into vapor or be boiled off 
by heat transfer. The time to pressurize is given by 

t = m.smg/Qt 

where: 

Qt is the volume flow rate based on the pressure drop, SCFM 

Table II gives the results of the saturation rule for a pressure drop of 20 
psig and pressurizing from 1 atm to 1.17 atm. 


Volume 

Mass, M 

Q 

T i me 

Gal Ion 

lbm " 

SCF 

Min 

100*000 

664 

1.27E5 

0.745 

200*000 

615 

1.18E5 

0.690 

300*000 

567 

1.09E5 

0.636 

400*000 

51? 

9.95E4 

0.582 

500*000 

471 

9.03E4 

0.528 

600*000 

422 

8.10E4 

0.473 

70 0.000 

374 

7.18E4 

0.419 

750*000 

350 

6.71E4 

0.392 


TABLE n Pressurization of the dewar from 1 atm at a pressure flow rate of 20 
psig. 
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Another contribution to the boiloff is the heat capacity of the dewar 
walls. On page 450 of Barron, an equation is given to determine the minimum 
thickness necessary for the walls of an elliptical shape. For a shpere, the minor 
diameter is the same as the major diameter. The equation for the thickness is: 

th = PDK/ (2S a E w - 0.2P) 
where: 


P is the maximum operating pressure 


D is the diameter, ft. 

K = 0.167 (2+ (D/Dl)l/2) and D/Dl is l for a sphere 
S a is the allowable stress, psi 
E w is the weld efficiency, per cent 


P is assumed to be 75 psia, D is 60.5 ft.. S a is 18,750, and E w is 100. Substituting 
these values into the equation, one obtains for the minimum thickness of 0.71 in. 
Using the density of steel of 7.84 g/cc and finding the net volume of the dewar wall, 
the mass of the inner wall is 9 . 3 x 108 g m . 

An equation for the internal volume of a sphere as a function of the 
height from the bottom of the sphere and the radius of a sphere is: 


V = 0.33Tth2 (3R-h) 


This equation can be written as a cubic equation in h with the radius as 30.25 ft. and 
V in gallons, the equation is 

h 3 - 90.75 h2 + 0.12766V = 0 


The surface area of the sphere at this height can* also be written as a function of the 
radius and the height. It is 
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S = 2rrRh 


Table HI shows the amount of hydrogen boiled off as a function of the 
volume of liquid initially in the sphere. The calculations iterate to find the height 
necessary to solve the cubic equation and then the surface area is calculated. The 
net difference in going from one liquid level to another in increments of 50,000 
gallons is given and the net surface area covered by liquid as it is filled to the next 
higher level is determined. The mass of steel knowing the wall thickness can then 
be calculated and the energy to cool it down 5 K is then arrived at. The mass of 
hydrogen boiled off can then be determined and the rate in SCFM is taken over 60 
minutes. 


V, Ui 

h, ft. 

A, f 12 

A, f 12 

n steel, 

a* Q, Joules 

*r H2 

SCFh 

50,000 

8.82 

1676 

1676 

2.32E. 

2.9E7 

3.25E4 

228 

100,000 

12*80 

2433 

757 

1.05E7 

1.31E7 

1 *47E4 

103 

150,000 

16*01 

3043 

610 

8.46E6 

1.06E7 

1.13E4 

Q3 

200,000 

18.84 

35ai 

538 

7.46E6 

9.33E6 

1 .04E4 

73 

250,000 

21.46 

<079 

498 

6 .9 1E6 

8 . o4E6 

9.67E3 

69 

300,000 

23.94 

4550 

471 

6.53E6 

8.16E6 

9.14E3 

64 

350,000 

26.34 

5006 

456 

6.33E6 

7 i 91E6 

3.36E3 

62 

<00,000 

28.68 

5453 

447 

6.21E6 

7.76E6 

8.60E3 

61 

<25,000 

29*85 

5673 

220 

3.05E6 

3.81E6 

4.27E3 

30 

<50,000 

31.02 

5893 

220 

3.05E6 

3.81E6 

4.27E3 

30 

500,000 

33*36 

6340 

447 

6*2 1 £6 

7.76E6 

8.60E3 

61 

550,000 

35.76 

6796 

456 

6.33E6 

7 . ?1E6 

8 . 86E3 

62 

L A A Art A 
u w V IVW 

TO A t 
JO. 4T 

7^ / / 
/ <iOO 

J V 4 
* 

6 • 53Eo 

5 . 16E6 

9 . 14E3 

64 

650,000 

40.86 

7765 

498 

6.91E6 

8.64E6 

? . 67E3 

68 

700,000 

43*69 

8303 

538 

7 . 46E6 

9.33E6 

1.04E4 

73 

750,000 

46*90 

8913 

610 

8 . 46E6 

1.06E7 

1.1SE4 

83 

800,000 

50*88 

9670 

757 

1.05E7 

1.31E7 

1 * 47 E 4 

103 

850,000 

60.25 

11346 

1676 

2.32E7 

2.90E7 

3.25E4 

228 


TABLE IQ . Boiloff from cooling down the heat capacity of the dewar. 
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The largest contribution to boiloff is from flashing of the higher pressure 
liquid to partial vapor into the low pressure dewar. An example of how to determine 
this is given in Van Wylen and Sonntag on page 122. (6) The First Law of 
Thermodynamics is used to find the state of the hydrogen in the dewar. The 
calculations were done twice, once using MKS units and once using English units. 
The results would be identical except that for the English calculation, I used the 
viscosity value of 9.54 x 10“6 lbm/ft-sec obtained from Scott. (4) In the MKS 
calculation, I used a value of 1.07 x 10~ 8 Pa-sec taken from a graph reprinted from 
Chapter II of Technology and Uses of Liquid Hydrogen . (7) The value 1.07 x 10“5 Pa- 
sec equals only 7.19 x 10 -8 lbm/ft-sec, and, therefore, the results do not agree 
exactly with each other. They are, however, in exact agreement when the friction 
factor is limited to a minimum value at high Reynold’s numbers 

The mass flow is given by: (8) 
m* = (ir2pAPD5/8fL)l/2 


Friction factors were needed in order to calculate the flow. Barron gives friction 
factor equations on page 135 for smooth pipes, but the equations he gives are only 
for Reynold's numbers less than 3 x 10 8 . (4) The reason for limiting the range for 
the validity of the equations can be seen from an examination of charts of friction 
factor versus Reynolds number as in "Flow of Fluids through Valves, Fittings, and 
Pipe, Technical Paper No. 410," (Crane Company), pp. A-24 and A-25. (9) 


For all pipes, the friction factor levels off to a constant value at high 
Reynold’s numbers. This constant minimum value depends on the relative roughness 
of the pipe. For very smooth pipes such as drawn tubing, the minimum will be 
0.0085 to 0.0095 for 4 inch tubing at Reynold's numbers of 3 x 10 8 and greater. But 
for commercial steel pipe, the minimum will be about 0.016 for Reynold's number 
greater than 2 x 10 8 . If the pipe in question has the roughness of commercial steel 
rather than drawn tubing, the Reynold's numbers of the mass flow are in this 
"minimum value" range and Barron's friction factor equations are not valid. 


The mass flows were calculated using several assumptions: 


A-12 



1) The friction factor equations are valid for all Reynold’s numbers, 

2) the friction factor reaches a minimum value of 0.0098 at Reynold’s numbers >3 
x 10®, 


3) The limiting minimum value is 0.016 instead of 0.0098. 

The analysis assumes that the hydrogen state point moves along the 
saturated-liquid curve as it is transferred, that is, the hydrogen is at the boiling 
point appropriate for the pressure at that location. We will assume that the tanker 
pressure is 22.5 psig and that the hydrogen enters the transfer line as a liquid (the 
quality or ratio of mass of vapor to the total mass flow is zero at the tanker end of 
the pipe). As the liquid flows along, the pressure falls, finally reaching 2.5 psig at 
the storage tank. Since the pressure falls, the boiling point decreases so some of the 
fluid becomes liquid at this new lower boiling point and lower enthalpy while part of 
it is vaporized to a gas, also at this new lower temperature but with an enthalpy 
greater by the heat of vaporization. There would be no net enthalpy change if there 
were not heat leak in the transfer line, no change in the potential energy of the gas, 
and no change in the velocity of fluid along its path. 


If there are elevations, some of the pressure drop goes into raising the 
potential energy of the fluid. We will show later that the heat leak of about 300W 
causes little additional vaporization of the hydrogen as compared to the pressure 
drop itself. Finally, since part of the fluid does change to vapor in the pipe, the 
mean density of the exiting hydrogen is less than that of the entering hydrogen. 
Thus, to conserve mass flow, the hydrogen must exit at a higher velocity than it 
enters. The increase in kinetic energy required ’’uses up" some of the pressure drop 
and leads to less liquid going to vapor. 

From Barron, 


(2 ( dh\ 

Jl IdP J 

X 2 = J - ±—L 


sat 


dP - V%-V? 

2 


- g (Z 2 -Z 1 ) + 


h %2 
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NOTE: I decided to use "V” for velocity and "v" for specific volume 


where: 

X 2 = quality of hydrogen at exit, 


dh 

dp sat 


change in enthalpy of the liquid for a change in pressure along the 
saturated-liquid curve, 


V 2 = velocity of exiting vapor and liquid, ft/sec. 

V 1 = velocity of entering liquid, ft/sec. 

(Z 2 -Z 1 ) = net change in elevation, ft. 

Q = heat leak into transfer line, BTU/hr. 
m = mass flow rate, lb/hr. 

hfg 2 = latent heat of vaporization for conditions at exit, BTU/lb 

In Figure 1, I have plotted (dh/dp) versus P for liquid hydrogen over the range of 
interest. The data are from NBS Monograph 94.(9). For this example, the transfer 
conditions will be assumed to be Pj = 22.5 psig = 2.53 atm, P 2 = 2.5 psig = 1.17 atm. 
Therefore AP = 20 psi = 1.361 atm. The first term in the numerator is just the area 
under the dh/dp curve over the region L17<P<2=53, This area was approximated by 
taking AP times the value of dh/dp at the midpoint, >^25.5 J/g atm, giving 34.7 J/g. 

We will assume net elevation changes to be zero. The velocity change 
can not yet be calcuated - we must iterate to get a final answer. Q= 1000 BTU/hr = 
293 Watts, so the last term in the numerator is 293 W/m. If nvr 8 kg/s, then Q = 
0.04 J/s and is completely negligible for this high mass flow rate. 

Then, as a first approximation, 
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X 2 = TTT^? = »-0783 


If, however, 7.83% of the transferred mass flashes to vapor, the mean density of the 
exiting liquid and vapor will be much less than the density of the liquid entering the 
transfer line. In order to conserve mass flow, the velocity of the exiting fluid must 
be greater than that of the entering liquid. 

p 2 = 1 = 1 

v m2 (I-X 2 ) vf 2 - X 2 v g 2 

where: 

v m2 = the mean specific volume at exit, 

Vf 2 = liquid specific volume at P = 2.5 psig, T = 1/kg 

Vg 2 = vapor specific volume at exit conditions, 1/kg 


Thus 

P2 = = — - — g/cm^ 

(0.9217) (14.254) + (0.0783) (652.9) 64.26 

P2 = 0.01556 g/cm^ = 15.56 g/1 = 15.56 kg/m^ 

Since the initial density was 66.21 kg/m3, the velocity change is larger. 

As v = m/(PiA), vi = 12.69 m/s for a mass flow rate of 8 kg/s. 

Thus V 2 = P1 vj = 4.255V j = 54 m/s 
P2 

2 2 

V2_ 2 Vi = 1377 J/kg = 1.377 J/g. 


The energy of 1.4 J/g is required to increase the momentum of the fluid 
and is supplied by part of the pressure drop. There is then less energy available to 
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provide the heat of vaporization of the liquid hydrogen. 


Our second approximation is kthen, 


X 2 = 


34.7 J/g - 1.4 J/g 
443 J/g 


0.0752 


We proceed in this manner and obtain the results shown in Table IV. We see the 
result converges to 0.0754 in three iterations. 


Iteration 

Kinetic Energy Change 
(J/q) 

(kcf/m 3 ) 

X 2 

0 

0.00 

66.21 

0.0783 

1 

1.38 

15.56 

0.0752 

2 

1 .29 

16.06 

0.0754 

3 

1 .29 

16.02 

0.0754 


TABLE IV: Results of iterative calculation for quality of hydrogen at exit. 

Discussion of flashing results 

The results show that the density of the exiting hydrogen is only 24% of 
the density of the entering fluid. The kinetic energy change required uses 4.85 psi 
of the 20 psi pressure drop. This implies that the mass flow calculations done in the 
first part of ther report could be greatly in error since, for a pressure drop of 20 psi, 
only '/'is psi are available to overcome the friction losses. But more importantly, 
the assumption that the flow could be calculated using only the properties of the 
fluid (density and viscosity) is invalid. In fact, whereas every gram of entering 
hydrogen consisted of 15.1 cc of liquid, every gram of exiting hydrogen consists of 
only 13.2 cc of liquid plus 49.2 cc of vapor. Moreover, the flashing results depend on 
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the validity of the homogeneous model which requires bubble flow. In the example 
in Barron the vapor was only about 20% of the exiting volume. In our case, it is 
nearly 80% of the exiting volume. The problem is really one of two-phase flow 
because of the large amount of flashing produced by the 20 psi pressure drop. 

Since the mass flow and flashing results are mutually interdependent, the 
results obtained above are not completely valid. The effect of flashing is to reduce 
the mass flow rate. The effect of the mass flow rate on the amount of flashing is to 
increase the amount due to the heat leak into the transfer line. Also, the kinetic 
energy term will be smaller, but as seen in Table IV, this is a small effect. Since the 
amount of flashing depends mostly on the total pressure drop (frictional pressure 
drop plus momentum pressure drop) , this value will probably not change much as 
long as the 20 psi drop is maintained during the transfer and as long as the flow 
remains mixed, i.e., the liquid and vapor don't stratify. 

If we use the average density, 41.12 kg/ m3, and assume the flow is still 
fully turbulent so that the friction factor has bottomed out at a value of 0.015, then 
after a few iterations we find 


P2 = 15.75 kg/m3, p aV g = 40.98 kg/m 3 
AP (frictional) = 17.7 psi, AP (momentum) = 2.3 psi 
m= 5.40 kg/s, and X 2 = 0.0771 

Therefore 416 g/s would be produced by flashing; corresponding to 


(416 g/sX60 s/m) 
2.362 g/ft 3 


10,570 SCFM 


where a standard cubic foot is taken to be at 1 atm and 70° F 

If the flow stratifies, the vapor and liquid may exit at different 
velocities. To handle that situation one might attempt to write down a set of 
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equations describing the flow at different segments along the length of the pipe and 
match the end points of each segment. In any case, the flowing hydrogen is going to 
undergo a pressure drop, so, for a 20 psig drop, ^7% - 8% of the mass of hydrogen 
leaving the tankers will arrive at the storage vessel as vapor. Knowing this, it is 
probably easier just to empirically determine the mass flow rate by measuring the 
time it takes to fill up the storage vessel. Then the rate at which gas will be vented 
during the transfer (as a result of flashing) can be calculated. 

Note that these results are a worst case since the tankers arrive with a 
pressure of only 5-10 psig and are then pressurized into 22.5 psig. The hydrogen is, 
therefore, not on the saturation curve when the transfer begins but is instead at 
some temperature below the boiling point at P =22.5 psig. To compute the flashing 
losses one needs to plot dh/dP versus P for the liquid in going from its bulk 
temperature in the tanker to its boiling point at 2.5 psig in the vessel. 

As the transfer proceeds, the liquid hydrogen in the tanker will warm 
through heat inleak, and the flashing losses will increase. If the heat transfer 
through the tanker insulation is known, the time required for the hydrogen to reach 
the boiling point can be calculated. (The heat capacity of the liquid hydrogen is 
about 11 J/ (g-k) for the conditions in the tanker.) 

The amount of liquid needed to cool down the transfer line depends upon 
what happens to the ’’hot" gas that exits the transfer line during this process. 
Perhaps there are valves which allow the gas to be vented directly until the transfer 
line is cold, and then the valves are reconfigured to allow the fluid to enter the 
storage vessel. 

If this is not the case, then one must compute how much liquid hydrogen 
will be boiled off when the storage vessel is not empty. If it is empty, then the 
amount of hydrogen required to cool the storage vessel down must be computed. In 
the former case, the amount of hydrogen that will be evaporated by the hot gas 
arriving at the beginning of the transfer will depend on the details of the 
construction. If there is a phase separator or vapor diffuser, the warm gas will be 
directed away from the surface of the cold liquid so that it vents before transferring 
much of its heat to the liquid hydrogen already in the vessel. 
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One of the last things needed to know in designing a hydride storage vessel is the 
temperature of the hydrogen entering the storage vessel. If the hydrogen is cold, 
this will provide some cooling to the storage beds. To find the temperature, the 
external heat gain from the ambient air must be calculated and equated to the 
heating of an equivalent mass of hydrogen. 

Three resistances are needed to find the external heat transfer rate. 
They are the external heat transfer coefficient, the internal heat transfer 
coefficient, and the resistance due to the conductivity of the pipe. The external 
heat transfer coefficient is assumed to be 20 BTU/hr-sq.ft-F for lack of better 
information. This is about midrange for a natural convection coefficient. The 
internal coefficient is more difficult to estimate. Since it depends on the 
temperature at which the properties of the hydrogen gas exist. If the bulk 
temperature of the hydrogen gas is used, then curves, iterations, or tables of values 
must be made to find where the external heat transfer from ambient temperature 
equals the sensible heating of the hydrogen. This varies directly with the mass flow 
rate of the hydrogen. 

The external heat transfer coefficient can be found from the following 
equation of the Nusselt Number: / 

N nu = 0.021 Np r 0 * 6 N re 0,8 
where: 

N nu Is the Nusselt number 
Np r is the Prandtl number 
and N re is the Reynold’s number 


The Nusselt number is equal to the heat transfer coefficient times the internal 
diameter of the pipe divided by the conductivity of the hydrogen gas. The Prandtl 
number for hydrogen is about 0.7 over a wide temperature range. Using this number 
and solving for the heat transfer coefficient: 
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h(t)= [0.01695 k(t)/D] [VD/v(t)/u(t)] 0.8 


where: 

h(t) is the external heat transfer coefficient 
D is the diameter, ft. 

k(t) is the conductivity of the hydrogen gas as a function of temperature 
V is the velocity of the hydrogen gas, ft/sec. 

v(t) is the specific volume of the hydrogen gas as a function of temperature, ft^/lbm 
u(t) is the viscosity of the hydrogen gas as a function of temperature, lbm/ft-sec. 

The external heat transfer is then equal to the following equation: 

Q = 1.8(300-t)/l/(2n0.434h(t)L)+l/(2n0.441(20L)+ln(1.016)/(2nk(t)L) 

Equations for the variables that are a function of temperature have been fitted to 
data using a least squares method. Since units for data varied, equations were fitted 
to the data as it was and then conversion factors were used in the final equation to 
get the correct result. 


Pt Ank er 

Pdrop 

fi*ss 

Mass 

Dens i t y 

Fr i c t i on 

Reynol ds 

Velocity 



F 1 ow 

F ) ow 


Fac tor 

Humber 


(pfti g > 

< ps * g > 

C 1 tom/ s > 

(kg/*) 

< ib »/ rt 3 > 



< f t /* ) 

13. 30 

13.00 

17.37 

7.97 

4.212 

0.0089 

6,490678 

40.710 

16.50 

14.00 

18.25 

8.28 

4. 201 

0.0089 

6742243 

42.396 

17.30 

13.00 

18.90 

8.37 

4.190 

0.0088 

698441 1 

44.832 

18.30 

16.00 

19.33 

8.86 

a m i 79 

A AAfifi 

7218048 

45 . $22 

19.50 

17. 00 

20 . 15 

9 . 14 

4.169 

0.0088 

7443896 

47.178 

20 . 30 

is. se 

20.74 

9.41 

4.158 

0.9987 

7662390 

48.692 

21 . 58 

19. 00 

21.31 

9.67 

4 . 1 ** 

0.0087 

7874680 

50. 139 

22.30 

20.00 

21.87 

9.92 

4.136 

0.0087 

8080643 

51 . 603 

23.50 

21.00 

22.41 

10 . 17 

4.126 

C . 0087 

8280907 

33.021 

24.58 

22 . 00 

22.94 

10 . 40 

4.115 

0. 0086 

8475838 

54.411 

25.50 

23.00 

23.45 

10.64 

4 . 104 

0.0086 

8665770 

35.776 

26.50 

24.00 

23.95 

10.87 

4.893 

0. 8086 

8851000 

57 . 118 

27.30 

23.00 

24.44 

11.09 

4.083 

0.0086 

9031794 

38.438 

28.30 

26.00 

24.92 

11.30 

4.072 

0. 0086 

9208394 

39.738 

29.30 

27.00 

23.39 

11.32 

4.061 

0. 0083 

9381018 

61.819 

39.30 

28.00 

23.84 

11.72 

4.050 

0.0083 

9 J 549865 

62.282 


1 

TABLE V. Result of calculation using English units. Friction factor fit assumed 
valid for all Reynold's numbers 
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Pt er 

Pdrop 

Mas s 

Mass 

liens i t y 

Fr i c t i on 

Reyno 1 ds 

Ve 1 oc i t y 



F 1 ow 

Flow 


F ac t or 

Number 


< ps i g ) 

<ps f g> 

Otg/») 

< 1 bm/ s ) 

( k g/m3 ) 



( m/s ) 

14.30 

12.00 

7.781 

17.15 

67.637 

0.0086 

8409846 

12.086 

15.30 

13.00 

8. 107 

17.87 

67. 463 

0.0086 

8762248 

12.623 

16.30 

14.00 

8.420 

18.36 

67. 293 

0.0886 

9100362 

13.146 

17.30 

15. 00 

8.721 

19.23 

67. 121 

0.0085 

9426193 

13.631 

18.30 

16.00 

9.01 1 

19.87 

66.948 

0.0083 

9740321 

14.142 

19.30 

17.00 

9.292 

20.49 

66. 776 

0.0085 

10,043,943 

14.621 

20.30 

18.00 

9.364 

21.09 

66.604 

0.8084 

10337924 

13. 088 

21.50 

19.00 

9.828 

21.67 

66. 432 

0. 0084 

10622998 

13. 344 

22. 30 

20.00 

10.084 

22.23 

66.239 

8.0084 

10899817 

13.990 

23.30 

21 . 00 

10.333 

22.78 

66.087 

0.0084 

1 1 168948 

16.428 

24.30 

22.00 

10.376 

23.31 

63.915 

0.0084 

1 1430895 

16.857 

25.30 

23.00 

10.812 

23.84 

63. 743 

0.0083 

1 1686107 

17.279 

26.39 

24.00 

1 1 . 842 

24.34 

63.370 

0.0083 

11934983 

17.693 

27.30 

23.00 

11.267 

24.84 

63.398 

0.0083 

12177884 

18. 101 

28.30 

26. 00 

1 1 . 486 

23.32 

63.226 

0.0083 

12413134 

18.302 

29.30 

27.00 

11.701 

23.80 

63.034 

0.0083 

12647030 

18.897 

30. 30 

28.00 

11.910 

26.26 

64.881 

0.0083 

12873838 

19.287 

31 . 30 

29.00 

12.116 

26.71 

64. 709 

0.0082 

13093803 

19.672 

32.30 

30.00 

12.317 

27.13 

64.337 

8.0082 

1 33 13133 

20.032 


TABLE VL Result when the friction factor is not limited even though fit may not be 
valid, ?4KS units. 


Pt Anker 

Pdrop 

Mass 

Mass 

Dens i t y 

Fr 1 c 1 1 on 

Peynel ds 

Velocity 



Flow 

Flow 


F AC tor 

Number 


<ps 1 g> 

<psig> 

< 1 Dm/s > 

<kg/* ) 

< 1 bo/ft 3 > 



<ft/*> 

13.30 

13.00 

16.72 

7.38 

4.212 

0.0098 

6,1 77996 

* 38. 746 

16.58 

14.00 

17.33 

7.86 

4.201 

0.0098 

6403020 

40.260 

17.30 

15.00 

17.91 

8. 12 

4. 190 

0.0098 

6619267 

41.727 

18.30 

16.00 

18.48 

8. 38 

4.179 

0.0098 

6827372 

43.151 

19.30 

17.00 

19.02 

8.63 

4.169 

0.0098 

7828641 

44.336 

20.50 

18.00 

19.33 

8. 87 

4.158 

0.0096 

7223078 

43.886 

21.30 

19.00 

20.06 

9. 18 

4. 147 

0.0096 

741 1403 

47. 203 

22.50 

20.00 

20. 53 

9. 32 

4.136 

0.0098 

7394877 

48.494 

23.30 

21.00 

21.03 

9.54 

4. 126 

0.0098 

7771493 

49.736 

24.50 

22.00 

21.30 

9.73 

4.113 

0.0098 

7944004 

58.994 

23.30 

23.00 

21. 95 

9. 96 

4. 104 

0.0098 

8111923 

32.208 

26.30 

24.00 

22. 39 

10.16 

•4.093 

0.0098 

3273331 

33.401 

27.30 

23.00 

22.83 

18. 33 

4.083 

0.0098 

8433077 

34.374 

28.50 

26.00 

23. 23 

18.34 

4.072 

0.0098 

8390789 

55.728 

29.30 

27.00 

23.66 

10.73 

4.061 

0. 0098 

8742871 

36.863 

30.50 

28.00 

24.06 

10.91 

4.030 

0.0098 

8,89 151 0 

37.983 


Table VII: Friction factor limited to values of 0.0098 or larger. In English units. 
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Pt ank er 

Pdr op 


Hass 

Bens i t y 

Friction 

Rtyno 1 ds 

Velocity 



Flow 

FI ow 


F ac tor 

Humber 


<ps * g> 

<ps i g> 

<kg^s> 

< 1 bm/s) 

(kg/».3) 




14.50 

12. 00 

7.295 

16. 88 

67.637 

0.0098 

7,8 8 3,6 2 2 

11.332 

15.50 

13.00 

7. 583 

16. 72 

67.465 

0.0098 

8197158 

11.810 

16.50 

14.00 

7.860 

17. 33 

67.293 

0.0098 

8495727 

12. 271 

17.50 

15.00 

8. 125 

17.91 

67. 121 

0. 0098 

8782651 

12. 718 

18.50 

16. 00 

8.381 

18. 48 

66.948 

0.0098 

9059037 

13. 152 

19.50 

17.00 

8.627 

19. 82 

66.776 

0.0098 

9325821 

13.575 

20.50 

18.00 

8. 866 

19. 55 

66.604 

0.0098 

9583806 

13. 986 

21.50 

19. 00 

9. 097 

20. 06 

66. 432 

0.0098 

9833685 

14.388 

22.50 

20. 80 

9, 322 

20.55 

66. 259 

0.0098 

10,076,059 

14.781 

23.50 

21 . 00 

9.539 

21.03 

66.087 

0.0098 

1031 1460 

15. 166 

24.50 

22. 00 

9. 751 

21.50 

65.915 

0.0098 

10540353 

15. 543 

25.50 

23. 00 

9.957 

21 . 95 

65.743 

0.0098 

10763154 

15.913 

26.50 

24.00 

10. 158 

22. 39 

65.570 

0.0098 

10980233 

16.277 

27.50 

25.00 

10.354 

22.83 

65.398 

0.0098 

1 1191925 

16. 634 

28. 50 

26.00 

10.545 

23.25 

65.226 

0.0098 

11398528 

16.986 

29.50 

27.00 

10, 732 

23. 66 

65.054 

0.0098 

11608315 

17.332 

30. 50 

28. 00 

10.914 

24. 06 

64.881 

0.0098 

11797533 

17.674 

31.50 

29.00 

11.093 

24.45 

64.709 

0.0098 

1 1990408 

18.011 

32. 50 

30.00 

11.267 

24.04 

64.537 

0.0098 

12179,145 

18.343 


TABLE VIP: Result when the friction factor is limited to values of 0.0098 or 
greater, in TviKS units. 


Pt Ank er 

Pdrop 

Hass 

Hass 

Dens 1 1 y 

Fr i c t 1 on 

Rtyno 1 ds 

Velocity 



Flow 

Flow 


F ac tor 

Number 


Cps 1 g> 

<ps i g> 

< Ibm^s) 

<kgxs> 

( 1 bm^f 1 3 ) 



<ft/s> 

15.50 

13.00 

13. 10 

5.94 

4.212 

8.0160 

484^323 

38.369 

16.50 

14.00 

13.58 

6. 16 

4.201 

0.0160 

5018697 

31.556 

17.50 

15.00 

14.04 

6.37 

4.190 

0.0160 

5188193 

32.706 

10*50 

16.00 

14.48 

6.57 

4. 179 

0.0160 

5351463 

33*822 

19.50 

17.00 

14.91 

6.76 

4. 169 

0.0160 

5509061 

34.907 

20.50 

18.00 

15.32 

6.95 

4. 158 

0.0168 

5661461 

35.966 

21.50 

19.00 

15.72 

7. 13 

4.147 

0.0160 

5809072 

36.999 

22.50 

20.00 

16. 1 1 

7.31 

4. 136 

0.0160 

5952251 

38.018 

23.50 

21.00 

16. 48 

7.48 

4. 126 

0.0160 

6091309 

38.999 

24 . 58 

22.88 

16.85 

7. 64 

4.115 

8 . 8 i 60 

6226523 

39.969 

25.50 

23.00 

17.21 

7. 80 

4.104 

0.0160 

6358139 

40.921 

26.58 

24.00 

17.55 

7. 96 

4.093 

0.0160 

6486375 

41.856 

27.50 

25.00 

17.89 

9. 12 

4.083 

0.0160 

6611428 

42.775 

28.50 

26.00 

18.22 

8.27 

4.072 

0.0160 

6733475 

43.680 

29.50 

27.00 

18.54 

8.41 

4.061 

0.0160 

6852677 

44,571 

30.50 

28.00 

18.86 

8. 55 

4.050 

0.0160 

6969,188 

45.449 


TABLE IX: Friction factor limited to values of 0.016 or larger in English Units. 
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Pt mnk er 

Pdrop 



Dens i t y 

Fr i c t ion 

Reyno 1 ds 

Velocity 



Plow 

F 1 ou 


F*c t or 

Humber 


c pft i g) 

( ps i g> 

> 


<kgs»3 ) 



(■/ft) 

i4.se 

12.ee 

5.718 

12.61 

67.637 

0.0160 

6,180764 

8. 882 

is.se 

13.ee 

3.944 

13.10 

67.463 

0.0160 

6424946 

9.237 

ie.$e 

i4.ee 

6. 168 

13.56 

67.293 

0.0160 

6638963 

9.618 

17.50 

is.ee 

6.368 

14.04 

67. 121 

0.0160 

6883856 

9. 969 

ie.se 

ie.ee 

6. 369 

14.48 

66.948 

0.0160 

7100488 

10.389 

ie.se 

17.00 

6. 762 

14.91 

66.776 

0.0160 

7309594 

10. 640 

2e.se 

ie. ee 

6.949 

15.32 

66.604 

0. 0160 

731 1804 

10.962 

21.se 

i9.ee 

7. 130 

13. 72 

66. 432 

0.0168 

7707659 

l 1 . 277 

22.50 

20. 00 

7,306 

16.11 

66.239 

0.0160 

7897632 

1 1 . 383 

23.50 

21 . 00 

7.477 

16.48 

66. 087 

0.0160 

8082139 

1 1 . 087 

24.50 

22.00 

7.643 

16.85 

63.915 

0.0160 

8261546 

12.183 

23.58 

23.00 

7.804 

17. 21 

63.743 

0.0160 

8436178 

\2 , 473 

26.50 

24.00 

7.962 

17.53 

63.370 

0.0160 

8606325 

12. 738 

27.50 

25.ee 

8. 1 IS 

17. 89 

63.398 

0.0160 

8772249 

13.038 

28.50 

26.00 

8.265 

18.22 

63.226 

0.0160 

8934183 

13.314 

29.50 

27.00 

8.411 

18.34 

63.034 

0.0160 

9092346 

13.383 

30. S0 

28.00 

8. 554 

18.86 

64.881 

0.0160 

9246926 

13.853 

31. 30 

29.08 

8.694 

19. 17 

64.709 

0.0160 

9398102 

14.117 

32. se 

33.00 

8.831 

19. 47 

64.337 

0.0160 

9546034 

14.377 


TABLE X: Result when the friction factor is limited to values of 0.016 or greater. 
MKS units. 
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J0 

* CaI c ul At e 

Mass F 1 ow 

20 

1 12 September 1983 

30 

OPTION BASE 1 


40 

INPUT ’’Which 

printer^" , Pr t 

30 

PRINTER IS Prt 

60 

D»4 . 334 * . 0234 


70 

L-700. *. 3048 


80 

Vi sc*. 1 *1 . 07E 

-4 

90 

Rho2* 1 . E*3*2 . 

01594/29.82 

100 

Rho3* 1 . E*3*2. 

01594/30. 98 

1 1 0 

fl«P 1 #D*D/4 . 


120 

PRINT USING ’’K"; “Pt4 nicer Pdrop 

ds 

Ve 1 oc i t y “ 


130 

PRINT USING “ 

K";" 

r " 
140 

PRINT USING “ 

K M ;“ (psi 9) ( psi 9) 


(m/s > " 



Masi Mass 

F 1 ow F 1 ou 

( k g/ s ) Cl bm^i > 


Density Friction 
Fac t or 

Ckg/m3> 


Reyno 1 
Nunbt 


130 PRINT 

160 IMAGE DDD . DD , 2X , DDD . D D f 2X , PD.DDD,3X,BD.DD t 2X,PDD.DDD,3X,Z. DDDD , 3X, 8D, 3X , DD 
D . DDD 

1?0 FOR Pc Ar* 1 4 . 3 TO 33 

180 Pdrope*Pc Ar-2 . 3 

180 Pdrop«6894 . 737*Pdrope 

200 Rho-Rho2*CRho3-Rho2>*< C < Pc Ar ♦ 1 4 . 69394 9 ) / 1 4 . 6 9394 8 > -2 > 

210 Re-1.E*5 

220 GOSUB Fric 

230 Mfl w*SQR<PI *PI ♦ Rho*D*D*D*D*D*Pdrop/'(8. *F*L> > 

240 REPEAT 

230 rtflu 0 «finw 

260 Re* 4. *Mf1u/<PI*V1 sc*D) 

270 COSUB Fric 

280 nn u*$QR<PI*PI *Rho*D*DeD*D*D*Pdrop/<8. *F#L> > 

290 UNTIL ABSCCMfl w-Mfl u0)/Mflw>< 1 . E-4 

300 Ve 1 oc i t y*Hf 1 w/ < Rho* A ) 

310 PRINT USING 160; Pc Ar ; P dr ope ; Mf 1 w ; MM u* 1 000 . '433. 3924 ; Rho ; F ; Re | Ve 1 oc 1 ty 

320 NEXT Pcap 
330 PAUSE 
340 Frici ! 

330 SELECT Re 

360 CASE <2300. 

370 GOSUB FI a« 

380 F«Fr_ 1 Am 

390 CASE 2300. TO 3000. 

400 COSUB FI Am 

410 GOSUB Fturb 

420 F«.3*(Fr lA*+Fr_turb) 

430 CASE 3000. ~T0 3 . E~6 

440 - - GOSUB Fturb 

430 F*Fr turb 

460 CASE ELSE 

470 BEEP 

483 GOSUB Fturb 

490 F*Fr_turb 

300 END SELECT 

310 RETURN 

520 FUsi ! 

330 Fr_1 Am«64 . /Re 
540 RETURN 
330 Fturbt \ 

360 Fr turb*. 0036*. 3*< (Re A <-. 32)) > 

370 RETURN 
580 END 


TABLE XI: Computer 
through X 


listing of program to calculate data presented in Tables V 

(MKS Units Version), 
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(jv<* p L 



P dh^dP 



(«tn) 


r^g)/atm 


1 

.1163 

32. 

189 



1 

.4230 

28. 

447 

-12. 

,2071 

1 

.8410 

23. 

461 

-07, 

, 1436 

2 

. 3406 

22 . 

, 57o 

-34. 

, 393 1 

2 

.9280 

21. 

136 

-83. 

.1203 

3 

.6135 

19. 

,7;3 

-32, 

,3723 

4 

.4033 

18. 

,630 

-01 , 

.3694 

3 

.3115 

17. 

,938 

-00, 

.7643 


Data. “Ploftc-ci ift r - t^. A. 

FIGDRE 1: Rate of change of enthalpy as a function of pressure differential. 
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Finally, the external heat transfer must be equal to the sensible heating 
of the hydrogen. The equation for the sensible Heat Q s of the hydrogen is: 

MgCp(t) (t— 30) 

where Mg is the mass flow rate of hydrogen 

Cp(t) is the heat capacity of the hydrogen gas as a function of temperature 

The following table summarizes the results. Most flows will be in the 
range of 1 to 1.5 lbm/sec. The temperature increase of the hydrogen for various 
lengths are given in Table XII. 


Mass 

Distance from 

Temperature 

Flow 

Dewar 

Rise 

Lbm/sec 

Ft. 

°K 

1.0 

100 

15.2 

1 .0 

200 

33.5 

1.0 

300 

51.2 

1.0 

400 

67.0 

1.0 

500 

79.5 

1.5 

100 

11.8 

1.5 

200 

24.5 

1.5 

300 

38.3 

1.5 

400 

52.0 

1.5 

500 

63.5 


TABLE XII: Temperature increase of vented hydrogen due to external heat 

transfer from the ambient air. 
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The total flow rate of hydrogen vented comes from the boiloff from 
heating up the sphere, ullage displacement, and flashing. Ullage displacement has 
not been calculated yet. It may be determined by the mass flow rate of liquid into 
the dewar and using the density of the liquid, determining the volume displacement 
rate. Using the density of the vapor leaving the dewar and the volume displacement 
rate, the mass flow rate of vapor leaving the dewar due to displacement can be 
calculated and an equivalent SCFM found. 


Sphere 

Sensible Heat 


Dewar L.FL Level 
Gal 1 ons 

Boi 1 -off 

SCFMx 1 0 J 
Ullage 

SCFMxl0 J 

Flashing 

Total Boil-off 
SCFMxl0 J 

50,000 

228 

1.903 

5.886 

7.783 

100,000 

103 

2.002 

6.196 

8.192 

150,000 

83 

2.056 

6.406 

8.456 

200,000 

73 

2.090 

6.706 

8.790 

250,000 

68 

2.115 

6.856 

8.965 

300,000 

64 

2.131 

7.136 

9.261 

350,000 

62 

2.145 

7.336 

9.475 

400,000 

61 

2.153 

7.456 

9.603 

450,000 

60 

2.163 

7.496 

9.653 

500,000 

61 

2.169 

7.636 

9.799 

550,000 

62 

2.185 

7.686 

9.855 

600,000 

64 

2.189 

7.716 

9.889 

650, U0O 

68 

2.133 

7.886 

10.063 

700,000 

73 

2.187 

7.896 

10.077 

750,000 

83 

2.191 

7.756 

9.941 

300,000 

103 

2.193 

7.606 

9.793 


TABLE XBI: Shows the contribution of each of the terms to the boiloff for a dewar 
at 2.5 psig and line pressure of 22.5 psig. 
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